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-p^ ■ Abstract 

^ In a general line with our first work [1], within hard thermal loop (HTL) approximation 
a general theory of the scattering for an arbitrary number of colorless plasmons off hard 
• thermal particles of hot QCD-medium is considered. Using generalized Tsytovich cor- 

^ ' respondence principle, a connection between matrix elements of the scattering processes 
^ ■ and a certain effective current, generating these processes is established. The iterative 

^ ' procedure of calculation of these matrix elements is defined, and a problem of their gauge- 
invariance is discussed. An application of developed theory to a problem of calculating 
energy losses of energetic color particle propagating through QCD-medium is considered. 
It is shown that for limiting value of the plasmon occupation number (~ l/g"^, where g 
is a strong coupling) energy losses caused by spontaneous scattering process of energetic 
particle off soft-gluon waves is of the same order in the coupling as other known losses 
type: collision and radiation ones. The Fokker-Planck equation, describing deccelera- 
tion (acceleration) and diffusion in momentum space of beam of energetic color particles 
scattering off soft excitations of quark-gluon plasma (QGP), is derived. 
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1 Introduction 



In the second part of our work we carry on with analysis of dynamics of boson excitations 
in hot QCD-medium at the soft momentum scale, started in jT] (to be reffered to as 
"Paper I" throughout this text) in the framework of HTL effective theory [2]. Here, we 
focuss our research on the studies of the scattering processes of soft-gluon plasma waves 
off hard particles of higher order within a real time formalism based on a Boltzmann type 
kinetic equation for soft modes. The nonlinear Landau damping process studied before 
jH] is simple example of this type of scattering processes. In a similar case of plasmon- 
plasmon scattering P], for sufficiently high energy level of the soft plasma excitations 
(exact estimation will be given in Section 5 below) all higher processes of plasmon - hard- 
partical scattering will give a contribution of the same order to the right-hand side of the 
Boltzmann equation. 

Our approach is based on the system of dynamical equations derived by Blaizot and 
lancu |1] complemented by famous Wong equation that describes the precession of the 
classical color charge Q = (Q"), a = 1, . . . , — 1 for hard particle in a field of incident 
soft-gluon plasma wave. Using the so-called Tsytovich correspondence principle intro- 
duced in Paper I with necessary minimal generalization for relevant problem, we establish 
a link between matrix elements of studied scattering processes and a certain effective 
currents, generating these processes. These effective currents appear in solving combined 
equation system of Blaizot-Iancu and Wong's equations in the form of expansion in powers 
of free gauge field A^^^"' and initial value of a color charge Qq. The coefficient functions in 
this expansion determines matrix elements of investigated scattering processes. The ap- 
pearance of classical soft-gluon loop corrections to the scattering processes, that in leading 
order in coupling can be formally presented in the form of tree diagramm with vertices 
and propagators in HTL or eikonal approximations, here, is a new interest ingredient. 

We apply the current approach to study of the propagation of high energy color parton 
(gluon or quark) through a hot QCD-medium and energy losses associated with this 
motion. Research of the energy losses of energetic partons in QGP at present is of a 
great interest with respect to jet quenching phenomenon jHl El E] and its related high-p^ 
azimuthal asymmetry [EllII], large-pt 7r° suppression [121113] etc., that have observed 
in ultrarelativistic heavy ion collisions at RHIC (see, review piP). Throughout the past 
twenty years by the efforts of many authors several possible mechanisms of energy losses 
have been analyzed: (l)elastic small distance collisional losses due to the final state 
interaction of high energy parton with medium constituent (Bjorken Braaten and 
Thoma [IE]); (2)polarization losses or losses caused by large distance collision^ (Thoma 

^Here, energy loss can be considered as a work performed by color charge polarizing of QCD-medium 
by own field. 
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and Gyulassy, Mrowczynski, Koike and Matsui J7], Braaten and Thoma ^Hl); (3) losses 
due to gluon bremsstrahlung induced by multiple scattering (Ryskin [TH]; Gyulassy, Wang 
and Pliimmer Baier at al. jlOl Ell ESj , Zakharov [SHj, Levin j^, Wiedemann [25] . 
Kovner and Wiedemann [221 etc.). The first two mechanisms often combined into one 
collision mechanism of energy loss, but it is convenient for our purpose to separated them. 

It was shown |B] that the elastic and polarization energy losses of partons in QCD 
plasma turned out to be too small for jet extinction, while induced radiative energy losses 
prove to be sufficient large to make itself evident in jet quenching in collisions of heavy 
nuclei^. For this reason recent theoretical stydies of parton energy loss have concentrated 
on gluon radiation [12] - [22] • 

However these 'traditional' approaches have difficulties accounting for the large jet 
losses reported at RHIC [^. This is the motivation for more detail analysis of mechanisms 
of energy losses already studied or their alternative formulations (in particulary, radiation 
losses) as well as the considered novel mechanisms. To the first case it can be assigned 
the works of Gyulassy, Levai and Vitev on construction an algebraic reaction operator 
formalism [2H], twist expansion approach developed by Wang et al. [21] • The work of 
Shuryak and Zahed [SD] considered synchrotron-like radiation in QCD by generalized 
Schwinger's treatment of quantum synchrotron radiation in QED, a work of E. Wang 
and X.-N. Wang [HI] connected with allowing for stimulated gluon emission and thermal 
absorption by the propagating parton in dense QGP, and also mechanism of coherent final 
state interaction proposed by Zakharov [32], can be assigned to the second case. 

In this work we would like to turn to study of just one more mechanism of loss (or 
gain) energy, the special and simplest case of which is polarisation loss. The approxima- 
tion within of which polarization loss is calculated, valids only for extremely low level of 
excitations for soft fields of medium. This is circumstance on which it was indicated by 
Mrowczynski in Ref. [2]- The expression for polarization part of energy loss derived in 
PTj up to color factors is exactly coincident with expression obtained early in the theory 
of usual plasma [SB] within standard linear response theory. However, we can expect that 
for ultrarelativistic heavy ion collisions generated QGP will be far from equilibrium in 
highly excited state. It can be valid at least for the subsystem of soft plasmons (Section 
2, Paper I), when interaction by soft waves preponderate over effects of particle collisions. 
One could study the influence of such off-equihbrium effects on parton propagation and 
radiation for unbounded QGP. For existence of intensity of soft radiation in medium, ad- 
ditional mechanism of deceleration (or acceleration) of energetic color particle connected 
with spontaneous and stimulated scattering processes of this particle off soft gluon ex- 

^As was mentioned in [^T] that although on theoretical estimations of the radiative energy loss of a 
hard parton much large than the elastic energy loss, a direct experimental verification of this phenomenon 
remains an open problem. 
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citations arises. In a limiting case of a strong gauge field |y4^(X)| ~ T, where T is a 
temperature of a system, this type of energy losses becomes of the same order in g as 
above mechanisms of energy losses and therefore can give quite appreciable contribution 
to total losses balance (this problem is discussed in more detail in Section 7 and 8 of 
present work). It is precisely here non-Abelian character of interaction of hard color par- 
ticle with soft-gluon field in QGP is to be completely manifested. The first work, where 
an attempt was made to accounting for interaction of energetic (massive) colored particle 
with stochastic background chromoelectric field, using the semiclassical equations of mo- 
tion, is the work of Leonidov [HI] . In this work an approach developed to the problem of 
stochastic deceleration and acceleration of cosmic rays in usual plasma has been used. 
However here, no account has been taken the fact that in the case of dense medium, what 
is QGP, in the scattering process, plasma surrounding traveling color charge not remains 
indifferent towards. In QGP the nonlinear polarization currents arise, essentially varying 
a physical picture. In the case of Abelian plasma this was shown by Gailitis and Tsytovich 
in . To account for polarization currents it is necessary to use a kinetic equation, that 
were not made in above-mentioned work by Leonidov. 

A more close to the subject of our research on a conceptional plane is the work of E. 
Wang and X.-N. Wang [3T] have been already cited above. Here, at the first a problem 
of influence on energy losses of stimulated gluon emission and thermal absorption by 
the propagating parton because of the presence of hard thermal gluons in the hot QCD- 
medium, was posed. This mechanism of energy losses is important in a medium with large 
initial gluon density (such it is proportional to gluon density) that it is really possible 
takes place by virtue of a strong suppression of high transverse momentum hadron spectra 
observed by experiments at RHIC ^3]. However, in Ref jSIj only gluon emission and 
absorption of hard thermal gluons (whose energy is of order T) with the use of the methods 
of perturbative QCD, is considered. For large density of soft gluon radiation consideration 
of contribution to energy loss of soft emission and absorption by the propagating parton, 
also becomes important. These processes (proportional to soft-gluon number density) 
by virtue of large occupation number of soft excitations adequately describe by using 
quasiclassical methods based on HTL approaches. 

Returning to the problem stated in this work, it is necessary run ahead note the follow- 
ing important circumstance: effective currents, determining the processes of spontaneous 
and stimulated scattering posses remarkable pecularity that they not explicitly dependent 
on mass of hard particle (in HTL-approximation), and therefore developed theory is suit- 
able equally for light color particle as well as massive ones. This in particular is reflected 
in that losses caused for instance by spontaneous scattering of energetic particle off soft 
excitations will be connected with not varying momentum of particle, but with rotation 
of its (classical) color charge in effective color space, governing by Wong equation. Here, 
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we have a principal distinction from energy loss of charge particles in usual plasma. In 
Abelian case as known 3^1, a rate of energy loss (or gain) for interaction with stochastic 
plasma field is inverse propotional to mass of hard particle, and therefore it is essential for 
light particles (electron, positron) and suppressed for heavy ones (proton, ion), although 
their polarization losses are practically identical. In non-Abelian plasma, at least within 
HTL-approximation the mass dependence can enter only through integration limits and 
a value of velocity of energetic parton, and therefore explicit suppression by mass of this 
energy loss type, is not arisen. In this connection note that study of propagation of heavy 
partons (c or b quarks) through QGP is of great independent interest. The first estimation 
of energy loss of heavy quark will be made by Svetitsky |37j for study of the diffusion 
process of charmed quark in QGP, and then the estimation will be given also in the works 
by Braaten and Thoma |16] Mustafa et al., Walton and Rafelski |38j (collision losses), and 
Mustafa et al. [SH], Dokshitzer and Kharzeev |40^ (radiation losses). The energy losses of 
heavy quarks are important not only with jet quenching phenomenon, but also for other 
important one for diagnostic of QGP: the modification of the high mass dimuon spectra 
from semileptonic B and D meson decays (Shuryak, Lin, Vogt and Wang jU], Kampfer 
et al., Lokhtin and Snigirev |l2]). 

Paper II is organized as follows. In Section 2 preliminary comments with regard to 
derivation of the Boltzmann equation taking into account the scattering processes of an 
arbitrary number of colorless plasmons off hard thermal particles, are given. Section 3 
presents a detailed consideration of the kinetic equation for the nonlinear Landau damping 
process, and here, also appropriate extension of the Tsytovich correspondence principle 
is presented. In Section 4 a complete algorithm of the succesive calculation of the certain 
effective currents, generating considered scattering processes, is given. In Section 5 on the 
based of the correspondence principle these effective currents are used for calculation of 
the proper matrix elements. Here, we estimate the typical value of plasmon occupation 
numbers, wherein one can restricted the consideration to acounting for only the contribu- 
tion from the nonlinear Landau damping process or all higher scattering processes should 
be considered. In Section 6 we discuss a problem of a gauge independence of matrix 
elements defined in previous Section. Section 7 is concerned with definition of general 
quasiclassical expressions for energy loss of energetic color particle for its scattering off 
soft plasma excitations. In Section 8 the energy loss caused by spontaneous scattering 
off colorless plasmons in lower order in powers of plasmon number density is analyzed in 
detail, and in Section 9 some problems connected with scattering of energetic particle by 
soft-gluon excitations lying off mass-shell, is discussed. Finally in Section 10 the Fokker- 
Planck equation, describing an evolution of a distribution function for a beam of energetic 
color partons scattering off plasmons, is derived. In the closing section we briefly discuss 
some interest and important moments concerned with mechanism of energy loss studied 
in this work, and remaining to be considered. 
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2 Preliminares 



As in Paper I we consider a pure gluon plasma with no quarks, where soft longitudinal 
excitations is propagated. Here, we restrict our cosideration to soft colorless exsitations, 
i.e. we assume that a localized number density of plasmons A'"'(p, x) = (A^p"^) is diagonal 
in color space 

where a, 6 = 1, . . . , A^ — 1 for SU{Nc) gauge group. We consider in Paper II the change 
of the number density of the colorless plasmons A^ as result of their scattering off hard 
thermal gluons. Wc expect in this case that the time-space evolution of scalar function 
A^p will be described by 

ON'' dN^ 

^ + • ^ = -Nlr^lNl] + (1 + Nl)r,[Nl] ^ -C[Nl], (2.1) 

where Vp = dcUp/dp is a group velocity of the longitudinal oscillations and uj^ = a;'(p) is 
the dispersion relation for plasmons. As it usually is, a functional dependence is denoted 
by argument of a function in square brackets. 

The more general expressions for generalized decay rate Fa and regenerating rate Fi 
can be written in the form of functional expansion in powers of the plasmon number 
density 

oo oo 

TdK] = E rf ^'^K], UK] = E rS'^^'^K], (2-2) 

n=0 n=0 

where 

^ = 1-^1 dT(2"+l) W2n+2{k\ P, Pi, . . . , P„; P„+l, . . . , P2n+l) (2-3) 

r(2n+l) J = /^/ dr(2"+^) W2nMM P, Pi, ■ ■ ■ , Pn; Pn+1, ■ ■ ■ , P2n+l) (2.4) 

X (1 + ATiJ . . . (1 + A^iJA^;.„,, . . . Nl^^^J^,[l + A]. 

Here, /k = /(k, x) is the disribution function of hard thermal gluons, and the phase-space 
integration is 



J dr(2'^+^) ^ / (27r)5(Ek + ^in - ^k+p,„-po.. - ^out) n (2-5) 



where 



^m = a;p+a;p^ + ...+a;p^, Pin = P + Pi 



^^out — ^p„+i + ■ ■ ■ + f^pan+l' Pout — Pn+1 + ■ ■ ■ + P2n+1, 

are total energies and momenta of incoming and outgoing external plasmon legs, respec- 
tively, and E]f_ = |k| for massless hard gluons. 

The (5-function in Eq. ()2.5j) expresses the energy conservation of the processes of stim- 
ulated emission and absorption of the plasmons. The function 

W2n+2 = W2n+2(k| p, Pi, . . . , p„; p„+i, . . . , P2n+l) (2.6) 

is a probability of absorption of n + 1 plasmons with frequencies u^, . . . ,uj'p^ (and 
appropriate wavevectors p, pi, . . . , p„) by a hard thermal gluon^ carrying momentum k 
with consequent radiation of n + 1 plasmons with frequencies cUp^^^ , . . . Wpj^+i (and the 
wavevectors Pn+i, • • • , P2n+i)- We note that for generalized decay and regenerated rates 
()2.Hj) and ()2.4|1 it was assumed that the scattering probability ()2.fj|l satisfies the symmetry 
relation over permutation of incoming and outgoing soft plasmon momenta 

W'2n+2(k| p, Pi, . . . , p„; p„+i, . . . , P2n+l) = W;2n+2(k| p„+l, . . . , P2n+i; P, Pi, • • • , P„). (2.7) 

This relation is a consequence of more general relation for exact scattering probability 
depending on initial and final values of momentum of hard particle, namely, 

W2„+2(k, k'l p, pi, . . . , p„; p„+i, . . . , P2n+i) = w;2„+2(k', k| p„+i, . . . , p2„,+i; p. Pi, . . . , p„). 

It expresses detailed balancing principle in scattering processes, and in this sence it is 
exact. The scattering probability ()2.(i|l is obtained by integrating of total probability over 
k' with regard to momentum conservation law: 

'W2n+2(k| p, pi, . . . , p„; p„+l, . . . , P2n+l) 2n5{E]^ + £in — -Ek+pin-pout ~ ^out) 
= j ^^'2n+2(k, k'l p, pi, . . . , p„; p„+i, . . . , P2n+l) 27r(5(Ek + Sin - ^k' - £oMt) 

dk' 

X (27r)35(k + Pin - k' - Pout) 

Such obtained scattering probability satisfies the relation ()2.7p in a limit of interest to 
us, i.e. when we neglect by (quantum) recoil of test particle. In the general case the 
expression ()2.7j) is replaced by more comphcated one (see Section 10). 

In Eqs. fj2.3p and ()2.4|) as in the case of pure plasmon-plasmon interaction (Paper I), 
we take into account the scattering processes only with equal number of the plasmons 

■^In the subsequent discussion a hard thermal particle of plasma, for which we consider the scattering 
processes of soft waves, we will call also a test particle. 
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prior to interaction and upon it, i.e. the scattering processes of the following "elastic" 
type 

g* + G ^ gt + G', for n = 0, 
g*+gt + G^g; + g* + G', forn = l, (2.8) 
. . . , 

where g*, gi, . . . are plasmon collective excitations and G, G' are excitations with charac- 
teristic momenta of order T. The scattering processes with odd number of the plasmons 
are kinematically forbidden by the conservation laws. Finally the scattering processes 
with unequal even number of incoming and outgoing soft external legs, i.e. the processes 
of "inelastic" type 

G^g*+gl + G', for n = 0, 

g* + gt + G ^ g; + g* + g: + g* + G', for n = 2, 

etc., have kinematic regions of momentum variables accessible by conservation laws not 
coincident with kinematic regions of the corresponding processes of elastic type (12. 8j) . and 
we suppose that a contribution of the last processes to the nonlinear plasmon dynamics 
of the order of our interest is not important. 

The scattering process for n = (Eq. ()2.8|) ) known as the process of nonlinear Landau 
damping in the case of a quark-gluon plasma was studied in detail in Ref. [3^ (we 
will consider this process in the following section in the somewhat different context). In 
the present work we would like to extend the approach developed in jS] to the scattering 
processes involving arbitrary number of plasmons. 

By using the fact that |k| ^ |p|, |pi|, . . . , |p2ri+i|, the energy conservation law can be 
represented in the form of the following "generalized" resonance condition 

^in - ^out - V ■ (pin - Pout) = 0, V = k/|k|. (2.9) 

In particular for n = we have resonance condition 

o;;,-^;,^-v(p-pi) = 0, (2.10) 

defining the nonlinear Landau damping process. Furthermore, one can approximate the 
distribution function of hard thermal gluons on the right-hand side of Eqs. ()2.3j) and ()2.4|1 

/k' ^ /k + (Pin - Pout) ■ (2.11) 

and set 1 + /k ^ 1 + /k' ^ 1 by virtue of /k, /k' ^ 1- 

We introduce the following assumption. We suppose that characteristic time for non- 
linear relaxation of the soft oscillations is a small quantity compared with the time of 
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relaxation of the distribution of hard gluons f^. In other words the intensity of soft 
plasma excitations are sufficiently small and they cannot essentially change such 'crude' 
equilibrium parameters of plasma as particle density, temperature and thermal energy. 
Therefore, we neglect by a space-time change of the distribution function /k, assuming 
that this function is specified and describe the global equilibrium state of non-Abelian 
plasma 

/k = 2 ^ 



Here, the coefficient 2 takes into account that the hard gluon has two helicity states. In 
the context of this assumption (as it will be further shown) the scattering probability 
'W^2n+2(k| p, pi, . . . P2n+i) dcpcuds upou the velocity v (a unit vector), but not upon the 
magnitude |k| of the hard momentum. This enables us somewhat simplify expression for 
collision term C[A^p]. 

We use the fact that the occupation numbers A^p. are more large than one, i.e. 1+Np_ ^ 
A^p.. Furthermore we present the integration measure as 



f dk 




fdQ^ 


1 (27r)3 ~ . 




Air 



where the solid integral is over the directions of unit vector v. Taking into account 
above-mentioned, considering ()2.11|) . the collision term can be approximated by following 
expression 

x{Sin - ^out)w'2n+2(v| p, Pi, . . . , p„; p„+i, . . . , p2„+l) A'pA'p^ . . . A'p^^^, • 

To derive the scattering probability 'W2n+2 in Section 5 we need somewhat different 
approximation of collision term. Setting 1 + N^, ~ N^, and /k' — /k, in the limit of a 
small intensity A^p — we have the following expression for collision term 

<^KlU'i-o==(/^|k|V.)|/f^^r.-«. (2.13, 



X W2n+2(v| p, Pi, . . . , p„; p„,+i, . . . , P2„+l)A'p^ • • • ^p, 



P2n + 1 



The kinetic equation ()2.1|) with collision term in the form of ()2.13|1 defines a change of 
plasmons number, caused by processes of spontaneous plasmon scattering off hard test 
gluon only. 
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3 Nonlinear Landau damping process 



In this section, we review the main features of the scattering probabihty for nonhnear 
Landau damping process derived in Ref. This will be done by another way (that has 
already used in Paper I) using Tsytovich correspondence principle ISH], admitting 
a direct extension to calculation of scattering probabilities for the processes of a higher 
order. We preliminary discuss a basic equation for a soft gauge field, that will play a 
main role in our subsequent research. We have already written out this equation in Paper 
I (Eq. (1.3. 8))"^. Here, it should be correspondingly extented to take into account the 
presence of a current caused by test particle passing through hot gluon plasma. 

One expects the word lines of the hard modes to obey classical trajectories in the man- 
ner of Wong r^] since their coupling to the soft modes is weak at very high temperature. 
Considering this circumstance, we add the color current of color point charge 

fQ'{x)=gv^Q^6^'\^-^rt) (3.1) 

to the right-hand side of the basic field equation. Here, is a color classical charge. 
Considering as a constant quantity we lead to the nonlinear integral equation for 
gauge potential instead of Eq. (1.3.8) 

*V-'^'^{p)At{p) = -J%',[A]{p)-fQ^{p), (3.2) 

where 

oo 

r/M{p) = Y.J^'^^'iA...,A), (3.3) 

s=2 

jl'^^A, • • • , ^) = ^ 9'-'J*^Z\::X(p^ -pi' • • • ' -ps)a''''Kpi)A'''''{p2) . . . A^^^-^ip,) 

s s 

x5^^\p-Y,Pi)'[[dpi, 

i=l 1=1 

and 

rQ'ip) = j^,v^Q'^6iv.p). (3.4) 

Here, we recall that the coefficient functions *r^^\ ''^''^ are usual HTL-amplitudes and 
*T>^^(j)) is a medium modified (retarded) gluon propagator in a temporal gauge defined 
by Eqs. (1.3.10) - (1.3.12). 

As in Paper I (Section 5) we consider a solution of the nonlinear integral equation 
()3.2j) by the approximation scheme method. Discarding the nonlinear terms in on the 
right-hand side, we obtain in the first approximation 

*V2_"'{p)Al{p) = -jj(p). 

■^References to formulas in ^ are prefixed by the roman number I. 
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The general solution of the last equation is 

A'.ip) = Af\p) - *VUp)jQiPl (3.5) 

where A^^^"'{p) is a solution of homogeneous equation (a free field), and the last term on 
the right-hand side represents a gauge field induced by test parton in medium. 

Furthermore, we keep the term quadratic in field on the right-hand side of Eq. ()3.2j) . 
Substituting derived solution ()3.5|) into the right-hand side, we obtain the following cor- 
rection to the interacting field 

4i)«(p) = -*:d^,(p)j(2)-(aw, A(°)) -*p^,(j9){j(2)-(A(°^ 

The last term on the right-hand side of this expression equals zero^ by virtue of f°''"^Q''Q'^ = 
0. The first term is associated with a pure plasmon-plasmon interaction and was analyzed 
previously in Paper I. Therefore now we have concentrated on new nontrivial terms putting 
in braces. Using explicit definition of the nonlinear current of the second order (Eq. ()3.3|) ) 
it can be shown that a second term in braces equals the first one. 

We determine thus a new effective current (more exactly, the first term in the expansion 
over free field) that is a correction to "starting" current ()3.4p . caused by interaction of a 
medium with color test particle 

4m"[^^°^](p) = 2J(2)'^^(AW, -*Pjq) (3.6) 

2 

= - ^ Q"^/t™ (p, -pi, -p + pi)A(°)'^^'^npi) "i^'^^Hp - piHi,siv ■ ip - p,))dp,. 

Here, in integrand the typical 5-function, determining the resonance condition ()2.10|) of 
the nonlinear Landau damping process arises. Using an explicit expression for effective 
current ()3.6|) one can define a scattering probability of plasmons off hard test particle. 
For this purpose according to Tsytovich correspondence principle it should be substituted 
Jq^"[A^°^] into expression defining the emitted radiant power of the longitudinal waves X' 
Eq. (1.4.4). However as in the case of the field equation (1.3.8) it needs to be preliminary 
carried out a minimal extension of an expression (1.4.4) taking into account a specific of 
considered problem. 

To the procedure of the ensemble average in Eq. (1.4.4) we add an integration over the 
colors with a measure 

dA 

f/Q= n^WQ"Q"-CA), dA = Nl-l, 

a=l 

^It is obvious that this term wiU be not equal to zero if the charges and Q'^ are refered to two 
different hard test particles. It will define contribution to effective current connected with the process of 
plasmon bremsstrahlung, that is a subject of research in our next paper j44j . 
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with the gluon Casimir Ca = N^. normahzed such that JdQ = 1, and thus 

JdQQ''Q' = ^6''\ (3.7) 

Besides it should be added an averaging over distribution of hard particles in thermal equi- 
librium: /(ik/(27r)'^/k . . . . Taking into account above-mentioned we will use a following 
expression for the emitted radiant power X' instead of (1.4.4) 

^---H.^(/*^V.)i,/.«/t^ 



X 



dpusign{uj)Q'^'^'ip) (J^;(v,p)J5^,(v,p))5(Re*A-i'(p)). 



Here, as distinct from (1.4.4) we remove an averaging over volume of a system, since 
in the definition of integration measure ()2.5|) a momentum conservation law is explicitly 
considered. Besides in notation ()3.8|) we take into account that for global equilibrium 
system an effective current (more exactly, the part caused by test particle) depends on 
momentum k only trough velocity v = k/|k|. According to corresponding principle 
(Paper I, Section 4) for definition of the scattering probability tt'2n+2(v| p, pi, . . . , P2n+i) 
it should be compared an expression obtained from ()3.8|) , with the expression determining 
the change of energy of the longitudinal excitations, caused by spontaneus processes of 
plasmon emission only 




^X) = (/SkVu) (3.9) 




£ /l7 / (^ / ^^'"""'^ ->-+2(v| p, p„ . . . , p„; p„^„ . . . , p,^^^)Nl . . . iVi^„_^^ . 



X 

n=0 ■ 



In derivation of the last equality the kinetic equation 1)2.11) with collision term in the limit 
of a small intensity A^p 0, (Eq. ()2.13|) ). is used. 

With all required formulas in hand now one can define the scattering probability for 
nonlinear Landau damping process. Substituting an effective current ()3.(ij) into ()3.8j) and 
comparing the obtained expression with the first term (n = 0) in the expansion on the 
right-hand side of ()3.9p . we obtain desired elastic scattering probability of plasmon off hard 
thermal particle. However such an obtained expression t(;2(v|p;pi) will not be a gauge 
invariant. This is associated with the fact that an expression for a first correction Jq^^ 
to "starting" current is not complete. As we have shown in Ref. [Sj, the effective current 
()3.6|) defines the scattering amplitude of a longitudinal wave off dressing the 'cloud' of a 
test particle. Diagrammatically this process of scattering is depicted in Fig.^ Here, the 
upper figure indicates what the Feynman diagram defines this scattering process - cutting 
of the effective self-energy graph before inserting a hard bubble along the gluon line. 
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I 




Figure 1: The process of the stimulated scattering soft boson excitation off test gluon through 
a resummed gluon propagator where a vertex of a three-soft-wave interaction is induced by 
*p(3)_ rpj^g blob stands for HTL resummation an the double line denotes hard test particle. 



There is a further process given a contribution of the same order as above- considered. 
This contribution on the classic language represents the normal Thomson scattering of 
a wave off thermal particle: a wave with the original frequency is a set in oscillatory 
particle motion, and an oscillating particle radiates a wave with modified frequency cjp^. 
Diagrammatically this process of scattering is depicted in Fig.|21 To define contribution 
of this scattering process to an effective current, account must be taken of influence of 
the soft random field on a state of color test particle itself. The particle motion in a wave 
field is described by the system of equations 

J^^,Qn)'i^. (3.10) 

Above, r is the proper time. The second equation is known Wong equation. Instead of 
the expression p.ip for the color current of a test charge now we must use more general 
expression 

The system of equations ()3.10|) is in a general case very complicated since it describes 
both the Abelian contribution to radiation connected with the change of trajectory and 
momentum of the test particle due to interaction with the soft fluctuating gauge field. 
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Figure 2: The Compton scattering of soft boson excitations off a test gluon. The upper diagram 
represents cutting of the effective 'tadpole' graph defining this scattering process. 

and the non-Abelian part of the radiation induced by precession of a color spin in the 
field of the incident wave. Our interest is only with a leading HTL-contribution over 
coupling constant in considered theory. This makes it possible to simplify the treatment 
and consider the hard test gluon as moving along the straight line, with constant velocity. 
In HTL-approximation the QED-like part of the Compton scattering is suppressed, and 
only the dominant specific non-Abelian contribution survives^. The processes of emission 
and absorption of plasmon is defined by a 'rotation' of color vector Q"^ = Q°'{t). Thus, in 
the limit of the accepted accuracy of calculation we need only a minimal extension of the 
color current of hard test particle 



^Let us note a similarity of this fact with cases which arize in considering of the problem of a different 
kind. In deriving energy loss of fast parton due to bremsstrahlung processes (radiation losses) it was 
shown that the QED-like part of the induced radiation associated with change of path, is suppressed in 
coupling constant as comparison with non-Abelian contribution, proportional to the commutator of color 
generators pil^ . 





(3.11) 



where a color charge satisfies the Wong equation 




(3.12) 
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Here, we turn to a coordinate time. 

To derive the oscillations of a color particle, excited by a soft random field and being 
linear in amplitude of field on the right-hand side of Eq. ()3.12|1 we set the color charge 
equal to its initial value: Q'^ito) = Qq, and we replace a field A''^{x) by a free field 

In this approximation the solution of the Wong equation ()3.12j) has the form 

Q'^[t) = Ql + igr'^Vo J ^ [e-*("-P)* - e-'^''-P^''']{v ■ A^''^\p))dp. (3.13) 

Substituting this expression into ()3.11|) and turn to Fourier-transformation we obtain 
additional to ()3.6p correction to "starting" current of test particle 

^■^'(P) = ' j£y3 r'^Qo I ^ [Siv ■ {p - Pi)) - 5{v ■ p)e-^(-^^)*°](t; • A^'^\p{))dp,. 

Here, the first term in square brackets in integrand defines a resonance condition for the 
process of the nonlinear Landau damping. The second term is associated with Cherenkov 
radiation that is absent in our case. In the subsequent discussion all terms in the expres- 
sions of ()3.13p type, containing initial time will be dropped, because they not contribute 
to the scattering processes of interest to us. 

Adding an obtained current jo^""^ {p) with current ()3.(ij) we derive complete expression 
for the first term in the expansion of an effective current in powers of a free field in the 
leading HTL-approximation 

J^qI\p) = j£y JkZ\Hp, -Pi)A^'^^^''^{pi)5{v ■ {p-p,))dp„ (3.14) 

where the coefficient function in integrand is defined as 

KZIMp, -Pi) = KZI'Mp, -pM = ^r^'Qo K,,, (v|p, -pi), (3-15) 
where in turn 

K^^, (v| p, -pi) = ^ + *V,,,,{p, -p,, -p + pi) r)'"'\p - p^)v,.. (3.16) 

The denominator v -pi in Eq. ()3.1(j|l is eikonal, that was expected in an approximation to 
the small-angle scattering of a high-energy particle. The sense of entering a symbol '(1)' 
over K^^^{y \ p, —pi) will be clear from the context below. 

We use the expression of an effective current ()3.14|1 for deriving of desired scattering 
probability W2(v| p, — pi). The procedure of a calculation of the scattering probabilities 
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with the use of Eqs. ()3.8|) and ()3.9|) will be given for general case in Section 6. Here, we 
present only the final result 



W2H P, Pi) = |T(v| p, -pi)|2 = \T^{p, -pi; 



(3.17) 



where the function 



T^ip,~Pi) = g' 



(1) 



i^MMi ~Pl] 



on— shell 



(3.18) 

represents the scattering amplitude of the nonlinear Landau damping process. The ex- 
pression ()3.17p was obtained in Ref. [3] by a different method, that will be mentioned 
below. 



4 The higher coefficient functions 

In previous Section we consider a simpler process of scattering of plasmons off hard 
thermal particle - the nonlinear Landau damping process. As in pure plasmon-plasmon 
interaction jjlj, a calculation of the scattering probability here, reduces to computation 
of some effective current (Eq. ()3.14|) ) generating this process. We have shown that this 
effective current apears in the solution of the nonlinear integral field equation ()3.2|1 . which 
defines interacting soft-gluon field in the form of an expansion in a free field A^f^ = 
^^(o)a^ and also initial value of a color charge Qq = (Qq). The last circumstance is a new 
feature of considered problem different from pure plasmon nonlinear dynamics. 

Based on results of Paper I and previous section we can rewrite now more general 
structure of an effective current in the form of a functional expansion in a free field A^^^ 
(and generally speaking, color charge Qo) generating the plasmon decay processes and 
the scattering processes of arbitrary number of plasmons'' by hard thermal particle 

jjtot)a[^(o)]^p) = Q^,v,6{v ■p) + Y: + E (4-1) 

l^^J s=l s=2 

where 

&A('^]{^,p) = ^ ^/ KZ::X(^\p, -pu . • • , -Ps)AW-^npi) ■ ■ • A('^^'^'{p,) 

s s 

x6{v- {p-J2P^))Il(^Pi^ 

i=l 4=1 

^Here, for concreteness we consider the scattering processes with longitudinal oscillation, but this refer 
also to transverse soft gluon excitations. 
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s\ 



fj.fll...fls 

s s 
i=l i=l 



The last sum on the right-hand side of Eq. ()4.H) represents an effective current generating 
plasmon decay processes. The current was studied in detail in Paper I, where a complete 
algorithm of succesive calculation of coefficient functions *r^^^^;;^''^, was proposed. Now 
our problem is in constraction of a similar algorithm for calculation of coefficient function 
^^'^l'.'.'^- general the structure of these functions is more complicated and tangled, as 
compared with *r^^^^' This follows from the fact that they represent infinite series in 
expansion over color charge Qq, i.e. 



(4.2) 



m=0 



where 



^mj^aai...as 

j^aai...as bbi...bm _ fifj.i...fj.a 



Qo=0 

Below we discuss a physical meaning of coefficients of this expansion on the example 
of particular computation. Here, we note that primary consideration will be focussed on 
computation of the terms linear in Qq in the expansion ()4.2j) that give leading contribution 
to the scattering processes of our interest. 

The more direct and explicit way of calculation of required coefficient functions K^^\= 
^^{'.'.'411) 'was presented in previous section. It is based on usual procedure of computing 
perturbative solutions of the classical Yang-Mills equation. However such a direct ap- 
proach for determination of an explicit form of the higher cofficient functions K^^\ s > 1, 
becomes very complicated and as consequence, ineffective. For deriving K^^'^ we use the 
approach suggested in Paper I, which was applied to obtain *r'^''^ with appropriate modi- 
fication. Remind that this approach is based on simple fact, namely, that the total color 
current '^'^ [A] = J^'^[yl] + jQ'[yl], entering into the right-hand side of the field equa- 
tion ()3.2|) . has two representation: by means of free and interacting fields, which must 
be equal each other. Thus in representation of free field A^^^ the current is defined by 
Eq. 1)4.111 . To rewrite explicit expression for total current in representation of interacting 
field we need to define the expansion for current jQ^[A]{p) similar to expansion ()H.H)1 . For 
this purpose we write the solution of Eq. ()H.12j) in the form 

Q''{t) = U^\t,to)Ql 

where 



U{t,to) = Texp{-ig f ■ A'^(r, vr))T'^dr} 

Jtn 
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-1 in '^tn J tn 



ft rri rTs-i 
dTi / dT2. . . 

'to J to J to 

is an evolution operator accounting for the color precession along the parton trajectory. 
Here, (T")'"^ = —if"-^'^, and in the last line we set v4^(t) = y4^(r, vr). Using this form of 
the solution of the Wong equation, we can present a current jQ^[A]{p) in the form 

fQl^](P) = Q>'^(^ ■ ^) + ^'^""'^'l^^t^ ^o) - 1]^"^^ ^. (4.3) 

Adding the obtained expression with nonlinear current ()3.3|1 . we derive a total expression 
for current in the representation of an interacting field 

OO 1 „ s s 

= E - 9^'' *r;:;v:X(p, -pu • • • , -p.)a"^'^Hpi)^"^''Hp2) • • • A'^^^^ips)5{p-Y.p,) n dp. 

s=2 i=l i=l 

(4.4) 

+ (2^ ^o"^^^^"^ ■ + "^^5 (2^ i (M^TT^ 

x(^; . A^^ip,)) ...{v A'^^{p.,))6{v ■ {p -j2p^))f[dp^ (T«^ . . .T^^'QI 

i=l i=l 

Here, in derivation of two last lines in Eq. ()4.4p we drop all terms containing initial time 
to- The interacting fields on the right-hand side of Eq. ()4.4|1 are defined by expansion 



A'^A'(p) = A(o)"'^(p) - *:D^^'(j9)ji*°*^"[A(o)](j9), (4.5) 

where the current j'^"^^"' is defined by Eq. ()4.H1 . Thus we have two different representations 
for total color current: Eqs. 1)4.111 and ()4.4p . which must be equal each other 

jJ*°*)"[A(°)](p) = jJ*°*)«[A](p). (4.6) 

Substitution of Eq. ()4.5p into the right-hand side of Eq. ()4.6p turns this equation into 
identity. As for pure plasmon-plasmon scattering P, for derivation of required coefficient 
functions K^^^ it should be differentiated left- and right-hand sides of equality 1)4.6)1 with 
respect to free field A^^^ considering Eq. ()4.4)) for differentiation on the right-hand side, and 
set = after all calculations. However in this case it is necessary to add differentiation 
with respect to initial color charge Qq to differentiation with respect to free field A^^\ 
and it should be added the condition Qo = to condition A^^^ = after all calculations. 
Below we shall give a few examples. 

The second differentiation of a total current yields 



^2j(tot)a[^](p) 



5Q^o M(o)«i/^i(pi 



5^^^i(v|p, -Pi] 



A(0)=0,Qo=0 



5{v{p-p,)) (4.7) 

?o=0 
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= {T'^^r' K^,, (v|p, -Pi) 5{v ■ ip-pi)), 

(1) iPT^ 

where Lorentz tensor Kf^f^^ is defined by Eq. (jH.lfjj) and represents a sum of two different 
contribution depicted in Figs.[T]and|21 

A more nontrivial example arises in calculation of the next derivative. It defines the 
process of nonlinear interaction of three wave with hard test particle 



53j(tot)«[A](p) 



A(o)=0,Qo=0 



S{v-{p-pi-p2)) 

Qo=0 

(4.8) 

[{T'^^T^^r' K,,,,, (v| p, -p,, -p,) + {T'^^T^-r^ K^,,^, (v| -p^, -Pi)] ^{vip-p.-p^)). 



5Ql (pi)M(o)<^2M2 



The color factor on the right-hand side of Eq. ()4.8|1 are multiplied by pure kinematical 
coefficients, which we will call partial coefficient functions, and are defined as follows 

-f^/.w/.2 -Pi,-P2) - 



{v-p2)iv ■ {P1+P2)) 

+ ^^^^^^^*^'''^'(P1+P2)T,,^,M2(P1+P2,-Pl,-P2) (4.9) 

V ■ {pi +P2) 

*T f,f,^^^^{p, -pi, -P2, -P + P1+ P2) *V'"''{p - Pi - P2)V^' 



~ , (1) 

+ *'^^,^,AP, -Pi, -P + Pi) {p - Pi) K^i^^ (v| p - pi, -P2) 

-*r^i.A(P, -P+P1+P2, -Pl~P2)*'D'"'' ip-Pl-P2)Vu*'D^^' iPl+P2)*Tyf,^^2iPl+P2, -Pl, -P2)- 

(1) 

Here, the fuction K ^^^^^ consists of five terms different in structure , whose diagram 
interpretation was presented on Fig. El Above-mentioned two examples suggest on that 
coefficient functions {5K1^^_^-;f^J 5Qq)qq=q have a color structure, similar to color structure 
of effective amplitudes *^1fff^',',°'ilf_f^, and thus coincide with color structure of usual HTL- 
amplitudes, first proposed by Braaten and Pisarski in Ref. j2] for an arbitrary number 
of external soft-gluon legs. To test this assumption we calculate a next derivative with 
respect to free field A^^^^: 



5Ql (pi)M(o)'^2M2(p2)M(o)"3M3(p3) 



(4.10) 

A(o)=0,Qo=0 

(1) 

^We note that with regard to an exphcit expression K (Eq. (|3.16(l ^ the last three terms in Eq. (|4.9|l 
can be also presented in another useful form 

*r;.piz.(p, -Pi, -P+Pi) *'D'"' {p~pi)v^, - *Tf,f,^^,^^{p, ~pi, -p2, -P + Pi +P2) {p-pi -P2)v^', 

where *f (4) is defined by Eq. (1.5.6). 
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Figure 3: The process of the sthnulated scattermg of thh'd soft boson excitations off hard test 
parton. 



Qo=0 



(1) 



[{T^-T'^'T^^r' ^MmM2/.3 (v| P, -Pi, -P2, -ps) + (perm. 1, 2, 3)] 6{v ■ (p - pi - P2 - Ps)), 



where 



(1) 
K 



(v|p, -Pl, -P2, -ps) 



VfjLi "^fJ-Z 



V ■ {Pl + P2 + P^) \{V-P2){V-{P2+P3)) 



+ v/V^XpI + P2 + p-s) *r^/^,^3^2(pi +P2 +P3, -Pl, -P3, -P2) 

+ v^*'D'"''{pi + Ps) *r^/^,^3 (pi + Ps, -pi, -P3)w^2 



1 



V^*V'"'\p2 + Ps) *rj,.^2^3(p2 + P3, -P2, -P3)f, 



Ml 



(1) 



W ■ (P2 +P3) 

- {*r^^.iz.(p, -Pl, -P + Pl) *I?'''''(p - Pl) i^i.>2^3 (v| P - Pl, -P2, -P3) 

+* r^^^A (P, -P3 , -P+P3 ) *^'''' {PsWu'^s ( V I I^Pl-P2 , -P3 ) *P (pi +P2 ) * F a'^^ ^2 (Pl +P2 , -Pl , -P2 ) 
+ *r^i.A(p, -p + Pl + P2 + P3, -Pl - P2 - Psy^'^'Xp - Pl - P2 - P3)Vu' 
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X*V^^'{pi +P2+ Psy^X'tii^it^-APl +P2+P3, -Pl, -P2, -Ps)} 
r ^ , (1) 

- |*r^/,,^2^(p, ~pi, ~P2, -p + pi+ P2) *'D'"' {p-pi- P2) K^,^^ (v| p-pi- P2, -ps) 

+ *r^A/.3i^(P> -Pl - P2, -P3, -P + P1+P2+ P3yi>'"''{p -P1-P2- P3)Vu' 

X*P^^'(pi +P2)*rv^,^,(pi +P2, -Pl, -P2) 
+ *^i,^L^Xu{p, -Pl, -P2 - P3, -P + P1+P2+ P3y'D'"''{p -P1-P2- P3)Vy' 
X*P^^'(P2 +P3)*rv^2/.3(P2 +P3, -P2, -P3)} 
- *r/,^i^2M3^(p, -Pl, -P2, -P3, -P + P1+P2+ P3) *V'"'\p -P1-P2- P3)Vy>. 

The right-hand side of this expression is represented for convenience as a sum of four 
groups of terms defined by derivation of currents jg^, J^^-*", J^^-*" and J^^"" (Eqs. ()3.2p . 
()3.3p ). respectively. Thus one can state that a color structure of the first term in the ex- 
pansion of arbitrary coefficient function ()4.2|1 linear over color charge entirely coincides 
with the color structure of usual (s + 2)-gluon HTL-amplitudes. 

In closing we consider the problem on physical meaning of higher power of color 
charge Qq in the expansion of coefficient functions K^^^^--""^^ (Eq. ()4.2p ). For this purpose 
we calculate the second derivative of the expression ()4.7|) with respect to Qq 



^3j(tot)a[^](p) 



(0^\3 

A(o)=0,Qo=0 



5{v ■ {p-pi)) 

?o=o 



9' ^j.,^j.,,ya, K,,,{y\p,-p,)5{v{p-p^)), 



where 



(2vr)= 

f ... (v| p, -Pl) = / U J'""'' ^ {v/V'^^UM 

J l(v-pi)(v- Pi + P 



- *P'^^'(pi + p'l) K,,^, (v| Pl + p;, -Pl) (4.11) 



V ■ (pi +p'i) 

+ *r^.i.2m (p, -K, -p + Pl + p'l, -Pl) (p'l)^;.,'*:^'^^'^^ (p - Pl - p'l)^^. 



(1) 



+ *r^,,., (p, -p;, -p + p'j ^^D^^'^i (p'i)t;.f ^P'^^'^^ (p - p'l) K,,^,^ (v| p - p;, -po j ■ p'jdp;. 

(2) 

It is easy to check by direct calculation that partial coefficient function K^^^^ (v| p, — pi) is 
symmetric relative to permutation p ^ — Pi, as it must be. The diagrammatic interpre- 
tation of different terms in the right-hand side of Eq. 1)4.111) is presented on Fig.EJ The 
presence of resonance factor 5(f ■ (p — Pi)) and integration over loop momentum p\ points 
to the fact that here, we are concerned with soft one-loop correction to nonlinear Landau 
damping process, that is suppressed by power of as compared with tree approximation 
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Figure 4: The soft one-loop corrections to nonlinear Landau damping process. 

()3.15|) . ()3.16|) . We note specially that the region of integration in loops is restricted by 
cone V ■ p'i = 0. The contribution 

y"*r^,^2^(p, -pi, -p + pi) *V'''''{p - pi) *VyXa{p - Pi, -Pi, -p + Pi+ p[) 

X *V^^'{p[)vx' "V^^Xp -Pi- p'lW ■ P'l) ■ {p - Pi))dp[ 

corresponds to the last term (in parenthesis) in Fig.|3J By virtue of a property of three- 
gluon HTL-amplitude: *Tf,^,^^.^{p, -pi, -pa) = -*r^^2;,i(p, -p2, -Pi), the integrand is odd 
under the interchange of p'l p — Pi — p'l and integrates to zero. 

It is clear that higher powers in Qo in the expansion ()4.2|1 is associated with soft loop 
corrections of higher orders, and they are suppressed relative to linear term (5ii'^^^/(5(5o)Qo=o 

(r) 

as ((7^)^ etc. The partial coefficients functions -ft'/xAti (v|p, —pi), r >2 (more precisely, an 
effective current connected with them) can be interpretated as "dressing" of initial cur- 
rent of hard test color particle or simple of a particle caused by interaction of this current 
with hot bath. For derivation of total probability of the process of the nonlinear Landau 
damping (process of elastic scattering for = in Eq. ()2.8p ). we must, generally speak- 
ing, summarize a series in expansion in a color charge degree Qo- It can be considered as 
a replacement of initial test particle by some effective quasiparticle on which in fact the 
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scattering process of soft-gluon wave takes place^. It is successfully that even in the case 
of highly excited gluon plasma (see next section) soft loop corrections are suppressed and 
therefore in leading order we can neglect by it. 

Finally we point to one interest relation, which exists between partial coefficient 

(2) (1) 

functions -ft'/x^ti — Pi), K ^^^^^ (v| p, — Pi, — P2) and four - soft-gluon effective subam- 
plitude *r^^j^2^3(p, — pi, — p2, — Ps) defined by equation (1.5.6). For its derivation we con- 
tract the expression (gH) with *P'^i'^i(pi)v, integrate over dpi and replace variables 
Pi — s> p'l, p2 — > pi and indices: /ii — z/, /X2 — /^i etc. After some algebraic transfor- 

(2) 

mations with use of Eq. (1.5.12) and explicit form Kfj_^-^ (v|p, — Pi) we obtain desired 
relation 

(2) f V V ~ I 

(v| P, -Pi) = - / 7 — , ,..2 {vy*V'"'{p + p\)v^,)5{v ■ p\)dp\ 
J [y ■ ypi + Pi) ) 

+ j K^,,, (v| p, -p;, -pi) *V^'\p^;)v,. b[v ■ p[)dp[ (4.12) 

+ ^J*i'^,u^,^x{p, -p'l, -Pi, -P + P1+ p'l) *V^^'{p -pi- p'i)vx' *'D'"'\p[)v^> 6{v ■ p[)dp[. 

The existence of such relation is not difficult for understanding by using diagramm rep- 
resentation of the functions entering in it. Thus the second term on the right-hand side 
of Eq. ()4.12|) corresponds to procedure of closing one of soft plasmon legs in Fig. El (in 
our case the legs with soft momentum pi) on hard test particle, as it depicted in Fig.HJ 
The last term on the right-hand side of ()4.12|1 is connected with closing two external soft 
legs of diagramm on the left hand side in Fig. I.l and attaching hard line of test particle 
to close loop. The physical meaning of the first term on the right-hand side of ()4.12|) is 
less clear (see however below). It is evident that relation of such a kind exists and for 
higher partical coefficient function and effective subamplitudes. It can be supposed that 
they are a consequence of gauge invariance (more precisely, covariance) of total effective 
current ()4.H1 with respect to gauge transformations concerned both a field A^^'>^{p) and 
a color charge Qq. We still return to analysis of relation ()4.12|) in our next work 03] in 
consideration of the process of nonlinear Landau damping with regard to rescattering of 
test particle off another hard test particle, where a reason of appearing the first term on 
the right-hand side of ()4.12|1 becomes more explicit. 

^Besides it can be interpretated classical correction loops as change of state not test color particle, but 
properties of a medium by the action of color field of test particle and interaction of initial "no dressed" 
current of particle with such modified medium. 
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5 Characteristic amplitudes of the soft gluon field 



In this Section we shall estimate for which typical amplitude of the soft gluon field, the 
contribution of the nonlinear Landau damping process to a generalized decay and regen- 
eration rates, rd[A^p] and ri[A^p], will be leading and for which value all terms in the 
expansions ()2.2j) will be of the same order in g. As a preliminary we define a general ex- 
pression for scattering probability W2n+2 = 'W^2n+2(v| p, pi, . . . , p^; p„+i, . . . , P2n+i)- Since 
the basic moments of computation of the scattering probability by a hard particle repeat 
reasoning used for deriving the scattering probability in pure plasmon-plasmon scattering 
in Paper I, Section 4, here, we restrict our consideration to brief scheme of calculation of 

W2n+2- 

Initial expression for deriving W2n+2 is Eq. (K18j) defining the emitted radiant power of 
the longitudinal waves XK The current entering into correlation function is given by 

n °° 

jQM^'^]iP) = (^^?0V(^-P) +E4?[^^°^]Kp), (5.1) 

where 

s s 

x6{v- {p-J2pi))]]_dpi. 

i=l i=l 

Here, in the coefficient functions K^'^^ we leave only the term leading order in g in the 
expansion ()4.2|1 . i.e. linear in color charge Qo 

KZ\:-:;M\p^ -pi, • • • , -Ps) ^ Kil\::l:\^\ p, • • • , -ps)qI 

In substitution of the current expansion (jS.ll) into a correlation function in integrand of 
Eq. ()H.8|1 we face the product of two series. As in deriving probability of plasmon decays, 
in this product it is necessary to leave only a sum of a product of terms having the same 
order in power of a potential A'^^\ In this case only a desired 5-function entering to 
integration measure / dT^'^"'~^^\ defining the generalized resonance condition ()2.9p . arise. 
Thus the emitted radiant power ()3.8p . taking into account emission caused by the processes 
of scattering plasmons off hard test particle, can be represented in the form of expansion 

oo 
s=l 

where 
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and 

^ - i (/^I^I'^a) l^Si i OSOS') (5.2) 

x/^ /dp ki i°'''Q'"''(p)A7™:--»(v| p, -pi, .... -p,)<;i::;f (viP, -p;, ■ ■ ■ , -pD 

s 



i=\ i=l 



The first term X*^°^ is connected with Cherenkov radiation (the linear Landau damping) 
which is kinematically forbidden in hot gluon plasma and therefore, this term can be 
setting zero. Let us consider the remaining terms X'-^'^\ s > 1. 

The integration over color charge in ()5.2|) is trivial. Furthermore we write out decou- 
pling of the 2sth-order correlator on the right-hand side of Eq. ()5.2|1 in terms of pair cor- 
relators by rule defined by us in Paper L Employing a condensed notion, Ai = A^^^'^^^pi) 
etc, we have 

{Af^ . . . Af->A^? . . . 4?)) = (slfiAf^A^^^) . . . (5.3) 

i=l 

^ (.!)^ n s^'^Q'^^'iP^) (- (^It^) {Ka^^ - o + Ni^A^^ + o}5(p. - p[). 

Here, in the last line, within the accepted accuracy, we replace equilibrium spectral den- 
sities Plpi) by off-equlibrium ones in the Wigner form (L3.15): P{pi) P{pi,Xi), slowly 
depending on x, futhermore, we take the functions /'(pi, Xj) in the form of the quasiparticle 
approximation (Eq. (L4.12)) and pass from functions /p. to the plasmon number density 

Substituting expression ()5.3|) into ()5.2|1 . performing an integration over 111=1 
taking into account the relation 

Siv-ip - ^p,))]'= ^r5iv{p- i^Pi))^ (5-4) 
i=i i=i 

we obtain instead of ()5.2|) 



i=l 
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27,5{v ■ {p - J2p,))T*'^^^-'^^\^r\ p, -pi, . . . , -p,)T»'^--=^(v| p, -pi, . . . , -p.) 



X 

1=1 



i=l 



Here, we introduce the function 

1/2 



«"-(v|p,-p„...,-p.)^,-fMP)) f^)x (5.6) 



X n ( ^ ) ( ^7== ) ^mT.:-.-;/ (-1 -Pi, • • • , -p.) Ion-Shell 



i=l \ ^"^Pi 



representing the interaction matrix element for the process of nonhnear interaction of s + 1 
soft longitudinal oscillations with hard test particle. Since external soft-gluon legs lie on 
the plasmon mass-shell, then account must be taken of the fact that resonance conditions 
()2.9p admit scattering processes with even number of plasmons, i.e. it is necessary to set 
in Eq. (IK311 s = 2n + 1, n = 0, 1, . . . . 

Furthermore, we multiply out terms in curly brackets in integrand of Eq. ()5.5|1 and 
use combinatorial transformation identical with transformation in Paper I. Confronting 
such an obtained expression for with corresponding terms in expansion one 

identifies the required probability W2n+2'- 

^^2n+2(v| p, pi, . . . , p„; p„+i, . . . , P2n+l) (5.7) 
1 (Pn) (Pn>Pi2) (piv,Pn) 

f rj-i*|a|fcrj-i|a|b I ^ ^ r-j-i*{a}br-j-i{a}fe _j_ ^ ^ r-j-i*{a}ferj-i{a}6 _j_ _j_ ^ ^ rjn*{a}6r-j-i{a}6 \ 



l<ji<n l<ii<i2<i 



where T^''^'^ = T'f"^^(v| p, pi, . . . , p„, — p„+i, . . . , — P2n+i)5 and for brevity we enter multi- 

(pii ) 

index notation {a} = (a, ai, . . . , a2„+i). The summing symbol denotes summing 

l<ii<n 

over all possible momentum interchange pi^, 1 <ii < n, by momentum (— pjj, n + 1 < 
ji < 2n + 1. The symbol Z]i<7i'<i2<n analogously denotes a summing over all possible 
interchange of momenta pair (pii,Pj2)5 1 < "^i < ^2 < by momenta pair (— pj^,— p^^), 
where n + 1 < ji<j2 < 2n + 1 etc. 

With explicit expression for scattering probability (Eq. ()5.7|) ) in hand, now we can turn 
to the estimation of typical amplitudes of the soft-gluon field. First of all we estimate an 
order of matrix element T""^ "2"+i''. In the soft region of the momentum scale the following 
estimation results from expression ()5.6j) 

rjnaai...a2n+lb ^ T^aa\...a2n+lh 
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The order of the coefficient function /^«"i - «2n+ife pg^j^ i^q estimated from arbitrary tree 
diagram with amputate 2n+2 soft external legs. Let us consider, for example, the diagram 
drawn in FiglHl The factor g is related to vertices, and eikonal propagator l/{v-p) ~ 1/gT 

1 2 2n+1 2n+2 



Figure 5: The typical tree-level Feynman diagram for scattering (2n + 2)-plasmons off hard test 
particle 

is related to an internal lines. From simple power counting of the diagram it follows an 
estimation 

T^aai...a2n+ib 2n+2 Ar{2n+l)/2 ^ ^1/2 ( r: o\ 

^MMi...M2n+i ~ 9 {gTy^+^ ^ ' ^ ' 

and thus 

rT^aai...a2n+ib ^ 2n+2 Ar(2n+l)/2 A/^ 

Using this estimation we derive from Eq. ()5.7p 

--+2~/"^^(^^iVf-^ (5.9) 
Furthermore integration measure has an estimation 

^^(2n+l) _ {^gTfn+2_ (5_^q) 

Power counting of the decay and regenerating rates ()2.3|) and ()2.4|) with regard to ()5.9|) 
and ()5.1()j) gives the following estimation 

■p(2n+l)aa' ^ j^2n+l rjn ^j^l ^n+l 

If we now set N^^ ~ — , p>0, then from the last expression it follows 

p(2n+l)aa' ^ g{2n+l) {2- p) ^2n+l rj. _ {^-H) 

For small value of oscillation amplitude (Eq. (1.6.5)) we have 

M2n+l)aa'\ _ ^2n+1^2n+ly n = 0, 1, . . . . 
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From this estimation it can be seen that each subsequent term in the functional expansions 
()2.2p is suppressed by more power of g'^. Thus for low excited state of plasma, correspond- 
ing to level of thermal fluctuations at soft scale jlHj , we can only restrict ourselves to first 
leading term in the expansions ()2.2j] 

describing the nonlinear Landau damping processes. 

In the other case of a strong field, |74^(X)| ~ T, for p = 2 from estimation (jS.llj) it 
follows 

(ri'^+')''''')^^^~iVf+iT, n = 0,l,... (5.13) 

i.e. the generalized rates are independent of g. All terms in the expansions ()2.2|) become of 
the same order in magnitude, and the problem of resummation of all relevant contributions 
arises. 

In closing this Section we compare estimations ()5.12j) and ()5.13j) with similar ones for 
generalized velocities in the case of pure plasmon-plasmon interactions. For low excited 
state of plasma, as shown in Paper I, the first leading terms in functional expansions of 
T^ i (according to Eq. (1.6.6)) have following estimation 

From this estimation and ()5.12|1 we see that four-plasmon decay process is suppresed by g. 
Such, when the soft-gluon fields are thermal fluctuations, the nonlinear Landau damping 
process is a basic process determining damping of plasma waves. 

A situation is qualitatively changed, when a system is highly excited. In this case we 
have an estimation for generalized rates, determining 2n + 2-plasmon decay processes, 
from Eq. (1.6.6) 

^p(2„+l)aa'^^^^^1^2ny^ ^ = 1,2,.... (5.14) 

Thus, from estimations ()5.13|) and ()5.14j) for sufficiently large intensity of plasma exci- 
tations, the processes of pure plasmon-plasmon interactions becomes dominant. There 
is certain intermediate level of excitations, when these two different scattering processes 
give the same contribution to damping of plasmon excitations. Considering the plasmos 
number density as A^p ~ 1/(7^ and comparing estimations ()5.11|) with (1.6.6), one can esti- 
mate roughly the value a, for which the orders of generalized velocities for these processes 
are comparable 

^ 3 1 (\YiN^\ 

^ = ^^2 + 2(1^]' ' 
this corresponds to |A^(X)| ~ (g/NcY^^T. 
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6 Gauge invariance of the matrix elements T^^^ 

Let us consider the problem of a gauge invariance of the interaction of matrix elements 
defined in previous section. From Eq. and expansions (jHZj), ()4.8|1 . ()4.1()|1 and ()4.11|1 
with respect to antisymmetric structural constant we see that a proof of gauge invariance 
is reduced to proof of the gauge invariance of functions presenting convolution of projector 
i'-fiip) = ■ " Pfj.)/{P ■ with partial coefficient functions taken on mass-shell^°, 

i.e. 

(r) (r) 

K{v\p, -pi, -ps) =Kf,^^,„^^ (v|p, -pi, -Ps)u^{p)u^'^{pi) . . .U^'{Ps)\on~-she\\■ 

{Q.l) 

We will show below that these functions are reduced to simple forms of (1.8.1) and 
(1.8.3) types, obtained in calculation of similar convolution with effective subamplitudes 

Let us consider the function 1)6.11) for s = 1 in linear approximation over color charge 

(1) ,rr™ 

Qo, i-e. for r = 1. Lorentz tensor K is defined by Eq. (jH.lfjj) . Using effective Ward 
identities for HTL-amplitudes and mass-shell condition, we derive 



K (v| p, -pi) = P^pI\— h *Tqq^{p, -pi, -p + pi) ^V""' {p - pi)v^\ 



• (6.2) 

on— shell 



It is easy to check that the term with a gauge parameter in Eq. ()6.2|) vanishes on mass- 

(1) 

shell. Now we consider a function K (v|p, — pi) in covariant gauge. For this purpose 
we perform the replacements (1.8.2) of projector and propagator on the left-hand side 
of Eq. ()fi.2)) . Then after analogous computations we lead to the same expression on the 
right-hand side of Eq. ()6.2p . Thus, we have shown that at least in the class of temporal 
and covariant gauges the matrix element for the nonlinear Landau damping process is 
gauge- invariant . 

We can assume that a similar reduction holds for an arbitrary partial coefficient func- 
tion. However, as in the case of a pure plasmon-plasmon interaction a proof of this 

(r) 

statement is imposible in general case by reason of absence of general expression -ft'^^i...^^. 
Similar to example in Paper I (Section 8) the only thing that we can make is to consider 
a contraction (jfj.lj) for three-point partial coefficient function, exact form of which ()4.9|l 

^"In Paper I we have used not entirely correct definition of a projector u^{p). It distincts from above- 
written by a sign. This in turn led to mistaken statement in Section 8, Paper I, on non-physical nature 
of matrix elements for plasmon decays with regard to odd number of plasmons. These matrix elements 
are related to actual physical processes, which would arise if they not kinematically forbidden by laws 
of conservation of energy and momentum by virtue of specific of spectrum of longitudinal oscillations in 
hot QCD plasma. 
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is known. Slightly cumbersome, but not complicated computations with the use of the 
effective Ward identities and mass-shell condition lead to the following expression 



(1) f 



V ■P2){V ■ {pi +P2)) 



+ —T^, r*roo.(-pi,-p2,Pi+P2)*:P'^'^'(pi+P2K' (6.3) 

V- {Pi+ P2) 

- *^mou{.P, -Pi, -P2, -p + Pi+ P2) *V'"'\p -pi- P2)v^i 

-. , (1) 

+ *root.(p, -Pi, -p^P\) {p - pi)K^,a{\\ p-pi,-p2) 

Toi.a(J5, -P+P1+P2, -Pl-P2Ti>'"'\p-Pl-P2)Vu*i>^^'{Pl+P2T'^X'00{Pl+P2, -Pl, -P2)j 



on— shell 



Here, it is also easy to check that all terms with a gauge parameter vanish on mass-shell. 
In deriving ()6.3|) we use relation p^K (v| p, — ]5i)|on-shoii = 0, that directly follows from 

If we now perform a replacements (1.8.2) on the left-hand side (l(i.3j) . then after anal- 
ogous computations we lead to the same expression on the right-hand side of ()(j.3|l . The 
function ()6.3p by virtue of Eqs. ()4.8|) and ()5.6|) defines the matrix element for scattering 
process including three soft-gluon waves and hard test particle. This process is kinemat- 
ically forbidden, i.e. the (5-function 

8{uj - uJi- UJ2- ^ ■ {V -Vi- P2)) 

has no support on the plasmon mass-shell. This suggest that all scattering processes 
involving odd number of plasmons are kinematically forbidden. In an expansion of the 
color current ()5.1|) all effective currents with odd number of free fields A'^^^ by 8- 
functions are equal to zero on mass-shell. The reflection of this fact is, in particular, a 
choice of the structure for generalized rates and Fi in collision term of kinetic equation 

(EH). 

As in Paper I one can note, that in spite of the fact that result ()(i.3j] is of only pure 
methodological meaning, nevertheless two examples ()6.2|) and ()6.3|) provide a reason to use 
considerably simple expressions of ()6.2j) type for all {2n + 2)-matrix elements ^°'°-i--°-2n+ib 
in particular calculations. 

All above-mentioned reasoning is concerned with a part of coefficient functions linear 
over color charge Qq. Genereally speaking, far from obviosly, that a reduction of ()(i.2|l and 
()6.3p types takes place and for partial coefficient functions relating to expansion terms 
()4.2j) for arbitrary order in Qq. Let us consider, for example, partial coefficient function 



4.1111 . that is connected with term in the expansion of -ft'^JJlj (v| p, —pi) quadratic over color 
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charge. Using effective Ward identities for HTL-ampIitudes and mass-shell conditions, we 
derive from ()4.11|) 



K,,, (v| p, -p,mp)u^^{p,) = pV, f {- \ {v,*'b^%[)v^, 



V ■pi){v ■ {pi +p'i)) 
^ - v,*V'"''{p, + p[) K,,o (v| pi + p[, -pi) (6.4) 



V ■ {pi+ p'l) 

+ *ro.,.,o(p, -P'l, -P + P1+ p[, -Pi) *V''^<ip[)v,,>*V''''''^ip - - p[)v,,, 
- To.,.,(p, -p[, -p + p[rV^<{p[)v,C'^''^'''^^p - p[) K,., (v| p - -pi)} ■ P'i)dp[. 

In deriving ()(i.4j] we use a replacement of integration variable: p'l —>■ p — pi — p'l and 
resonance condition v ■ {p — pi) = 0. It is not difficuh to check that all terms in integrand 
()6.4p with a gauge parameter vanish on mass-shell. Performing a replacement (1.8.2) on 
the left-hand side, after anologous computations we lead to the same expression on the 
right-hand side of Eq. ()6.4|1 . Above-mentioned examples ()6.2|1 . ()6.3|1 and ()6.4|1 suggest 

(r) 

that all terms -ft'^/ij, r > 2 in the expansion of -ft'^^j and higher coefficient functions 
j^aai...a2„+i ^ > Q causc gauge- invariant convolution in the form (16.11) . 

At the end of this section we would like to discuss a possible existence of very nontrivial 

connections between effective subamplitudes *r^^^...^4^^3, defined in Paper I and partial 

(1) 

coefficient functions linear over color charge K^^-^,,,^^^^^^. As was mentioned at the end 
of Section 3 the matrix element for nonlinear Landau damping process was obtained in 
[3 by an alternative method distinct from the method proposed in Sections 3 and 4 of 
present work. In it was shown that the obtained nonlinear Landau damping rate for 
longitudinal modes 

If \ O 2Ar fjhh.rl ( Z;(p) Z;(pi) ^ ~ 

7 (p) = - 2g N^j ^ N.\ ^au^^^^ 2^ u-{p,) ] ^ r(p,p„ -p, -p,), 

\ I \ \i \L I PI \j i / / on— shell 

where 

*f (p,pi, -P2, -Ps) = *V^^^^^^''{p,Px, -P2, -P3)Uf,{p)u^,{pi)Uf,2{P2)u^,{p3)\on-she\h 

and the function *r'^'^^'^'^'^^ is defined from Eq. (1.5.6) by direct transformation of imaginary 
part from convolution of *r(p, pi, —p2, —ps), can be introduce also in the form 



J (27r)3 P'\2ulu^{p)2ujiy{pi)J^^^_^^^ 
x(^p-^pj/^<^(^i,-^pi - V - (p-pi))|-ft^ (v|p, -pi)| 
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where the function 

(1) (1) 

K{v\p, -pi) = K^^^ (v| p, -pi)u^{p)u^^ (pi) Ion-shell 

and m^i is Debay screening mass. Thus an existence of two different representations of 

(1) 

nonhnear Landau damping rate is caused by relation between the functions *r and K: 

Im*f -p2, -p3)|p2=P,P3=Pi (6-5) 

= - K - - - V • (p - Pi)) I (v| p, -pi) |2. 

Based on obtained relation ()6.5|) we can make speculation relative to connection between 

(1) 

functions *r and K of arbitrary (even) order in external soft lines 

lm*T{p,pi, . . .,P2n+l, -P2n+2, • • • , -P4n+3) |p2„+2=P,P2n+3=Pl,-,P4„+3=P2„+l = (6-6) 
-nml,{£i^-£ont)j S{Sir,-Sout-V-iPm-Pout))\K {y\p,Pi, . . .,Pn, -Pn+1, -P2n+l))P- 

Here, 

T(p,pi, . . . ,P2n+l, -P2n+2, • • • , -Pin+s) 

= *f^A.i.../.4„+3(P>Pl> • • • ,PAn+3)u^{p)u^'{Pl) ■ • •M^'"+'(P4n+3) Ion-shell, 
(1) 

and the function K{-v\p,pi, . . . ,p2n+i) is defined by Eq. ()6.1|) . 

The equality is actually reflection of relation between the imaginary part of 

self-energy n(a;) and the interaction rate r(a;) in the form proposed by Weldon [46 j 

lmU{uj) = -uT{u). (6.7) 

In our case we have 

dp.Nl^ T(p, p„ -p, -p,) cuTiu;) ~ J dp,{u;l-ulJ J^\K (v| p, -p,) \\ . . . 

We remind that formula ()(j.7|l was obtained by Weldon from analysis of Feynman diagrams 
with the use of bare propagators and vertices. Therefore nontrivial moment in relation 
()6.5p is the fact that it involves resummed propagators and vertices. The example of 
relation of such a kind can be found in the work of Braaten and Thoma JH] ■ 
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7 Energy loss of energetic parton in HTL approxi- 
mation. Initial equation 



As an application of the theory developed in previous Sections we study a problem of 
calculating energy loss of energetic color parton^^ traversing the hot gluon plasma, i.e. 
energy loss due to a scattering process off soft boson excitations of medium in the frame- 
work quasiclassical (HTL) approximation. As initial expression for energy loss we will use 
a classical expression for parton energy loss per unit length being a minimal extension to 
a color freedom degree of standard formula for energy loss in ordinary plasma [HS] 

r/2 

-^ = ^Jim- j Jd^dtJdQReCrQ{^,t)-E^Q{^,t)) (7.1) 

X |v| r_^^^ 

= — Jim ^-^ J dpdu J dQ Re (Jq (p, uj) ■ E^(p, u)). 

Here, v is the velocity of energetic parton. Check from above is introduced to distinguish 
the velocity of energetic parton (color particle which is external with respect to medium) 
from the velocity v (|v| = 1) of hard thermal gluons. The sign on the left-hand side 
of ()7.H1 corresponds to the choice of a sign in front of the current in the Yang-Mills 
equation (1.3.1). Chromoelectric field Eq(x, t) is one responsible for parton at the site of 
its locating. In our accompanying paper |44| the expression ()7.H) will be extended also 
to a case of energy loss due to plasmon radiation induced by scattering in a medium off 
hard thermal particles (plasmon bremsstrahlung) . 

In the expression 1)7.11) as a color current Jq(x, t) one necessary takes an effective 
dressed current of energetic parton, that arises as a result of a screening action of all 
thermal hard particles and interactions with soft color field of plasma. As resulted from 
an expansion of coefficient function ()4.2|) . this effective current represents infinite series 
in expansion over color charge Qq of energetic parton 

oo 

Here, the term linear over color charge, i.e. the term with m = is leading in cou- 
pling constant. Let us write once again the general expression of the effective current 
jQ^[yl'^'^^](p) for the convenience of further references 

oo 
s=l 



"^As energetic color parton (or simple parton) we take either energetic quark or gluon. 
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where Jq]^^{p) = g / {2'kY5"'^ v^5{v ■ p), v = (l,v) is initial current of energetic parton, 
and 

x5{v■{p-j2P^))f[dp^. (7.3) 

i=l i=l 

The chromoelectric field caused by the current ()7.2|1 is defined by the field equation in the 
temporal gauge 

E-(p) = -^uJ*V^'^ip)J^i[A^'^]ip)Ql 
where the soft-gluon propagator in a given gauge reads 

*V^^ip) = (^^^ ^*^\p) + [s^' - ^)*A*(P), (7.4) 

*V'°{p) = *V'^\p) = *poo(p) = 0. 

Substituting this expression for field Eq{p) in Eq. ()7.H) and integrating over color charge 
by using ()3.8|) we lead to formula for energy loss, instead of ()7.ip 



^ = ^ Jim / dpdulm {r^fip) *V^^{p)Jt{p))- (7.5) 

Substituting the expansion ()7.2|1 into Eq. ()7.5|1 we will have in general case the terms 
nonlinear in soft fiuctuation field A^^^'^{p) (more exactly, their correlators). These non- 
linear terms define a mean change of parton energy connected with existence of both 
spatially-time correlations between fiuctuations of soft gauge field, and correlations be- 
tween fiuctuations of a direction of a color vector of energetic parton and fiuctuation of 
soft field of system. The existence of such correlations result in additional (apart from 
polarization P^) energy loss^^ of moving particle. 

First of all we write the expression for energy loss connected with initial current 
Jqj^^ip)- Substituting this current into ()7.5p and taking into account a structure of a 
propagator ()7.4|) we obtain 

(-^y =-^(^^yrfp^{lm(/*A'(p)) + (vxp)2lm(*A*(p))}5(t;-p), (7.6) 

where = g'^/4:7c. In derivation of ()7.6|) we use a rule ()5.4j) . This expression defines 
polarization loss of energetic parton, connected with large distance collisions. The con- 
clusions about energy loss in works ^Zj are valid, generally speaking, only for equilibrium 

^^This kind of loss sometimes is called fluctuation loss (see, e.g. 
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(nonturbulent) plasma. However for sufficiently high level of plasma excitations (strong 
turbulent plasma) contributions to energy loss connected with higher terms, Jq^"^, s > 1, 
in the expansion of effective current, become comparable with polarization loss ()7.6|1 and 
therefore these contributions are necessary for considering. This can be seen from the 
level of effective current ()7.2|) . Really, let us write an estimation of oscillation amplitude 
of soft field A^-^^"'{p) in the form 

By estimations (1.6.4) and (1.6.5) a value of parameter 6 = corresponds to the case of 
highly excited state of gluon plasma, and the value 5 = 1/2 corresponds to weakly excited 
state or level of thermal fluctuations at the soft momentum scale From expression 
for currents ()7.3p and estimation for coefficient functions ()5.8|) it is follows 



From these estimations it can be seen that for a small value of oscillation amplitude, 
i.e. for 6 = 1/2, each subsequent term in the expansion ()7.2j) is suppressed by more 
power of g^^"^. Here we can restrict ourselves to the ffist leading term Jq^^^ip) defining 
polarization loss ()7.6|1 . and higher order terms in expansion of the effective current will 
give perturbative corrections to expression ()7.6|) . In another limiting case of a strong field, 
when 5^0 from these estimations it follows that all terms in the expansion ()7.2|) become 
of the same order and correspondingly are comparable to value with polarization one. 

In what follows we will suppose that a value of parameter 5 is different from zero, but 
it is sufficiently small to consider first two terms after leading one in the expansion of the 
effective current 

Substituting the last expression into Eq. ()7.5|) we obtain the following after ()7.6p term in 
the expression for energy loss of energetic parton, that we write as a sum of two different 
in structure (and a physical meaning) addends 



Here, the ffist term on the right-hand side 

/ dE^^'^Y 1 Ca (27r)4 



dx 
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defines energy loss due to the process of spontaneous scattering^^ of energetic parton off 
plasma waves (i.e. plasma excitations lying on mass-shell). The second term 

+(i£""(p)-i""'(p)4°i°'(p)>}, 

as will be shown in Section 9, is different from zero for plasma excitations lying off mass- 
shell. The integrand in Eq. ()7.9|1 is proportional to 6{v ■ p), that gives ground to assign 
the second term to polarization loss ()7.fij) . more exactly, to its correction due to nonlinear 
effects of medium. Hereafter the expressions ()7.8|) and ()7.9|) for brevity will be called 
diagonal and off-diagonal contributions to the energy loss connected with diagonal and 
off-diagonal terms in a product of two series ()7.2|1 . 

Our main attention is concerned with an analysis of expression of energy loss due to the 
processes of spontaneous scattering off plasma waves, i.e. with expression ()7.8|) . assuming 
thereby that this contribution to energy loss is a main in general dynamics of energy losses 
in this approximation. By using the expression ()7.4|) of equilibrium soft-gluon propagator, 
the diagonal contribution to energy loss similar to ()7.(i|l reads 

+ {l(pxjS"°''(P))r)In'CA'(p))}. 

Following by common line of this work, the contribution to energy loss caused by 
spontaneous scattering off longitudinal plasma waves (plasmons) is of our interest, i.e. on 
the right-hand side of Eq. ()7.10|) we leave only contribution proportional to Im(*A'(p)). 
By using an explicit expression for current Jq]^^{p)-, Eq. ()7.H|1 and also a definition of a 
coefficient function Eq. the diagonal contribution can be introduced in the form: 

(-'-^J- (^) (4) <^i- 

x\m{*^\p))6{v-{p-pi)), (7.11) 

where 

K''' (v|p, -pi) = + *r''\p, -pi, -p + pi) "V"' {p - pi)v^' . (7.12) 
V ■ Pi 

^^The notion of spontaneous and also stimulated scattering will be considered in Section 10 in more 

detail. 
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In deriving ()7.1H) in the spectral density P^^{pi) we leave only a longitudinal part 
(Eq. (1.3.16)) 

UJl pi 

and going to an average value of the chromoelectric field squared 

/'(Pi)--A(Ef), 



UJ 



2 IVl- 
1 



The expression ()7.1H1 is general, since it takes into account the availability in the system 
of plasma fiuctuations lying both on plasmon mass-shell and off-shell. To define energy 
loss due to the scattering of hard parton off plasma waves in integrand on the right-hand 
side of Eq. f)7.11|) it should be set 

Im(*A'(p)) ~ -7rsign(cu) 5{Re*^'^\p)) = -7rsign(cu) [5{uj - J) + 5{uj + uj% 

Substituting these expressions into ()7.1H) . integrating over du and duJi after some algebraic 
transformations we lead ()7.1Hl to the follows form 



dx I IvI V 27r^ 



J dpdpiW2{y^\ p, pi) ulNl^Siv ■ {p - pi)) (7.13) 



dpdpiW2{ir\p,Pi){ujlNl^ + ul^Nl}6{v ■ {p - pi)). 



27r2 

Here, a scattering probability W2(v|p, pi) is defined by Eqs. fl3.17|l . ()3.18|1 . where the 
factor g^Nc enters to one on the right-hand side of Eq. ()7.13j) . and a replacement v ^ v 
is made. In deriving ()7.13p we drop a term with probability W2(y \ p, — Pi), that defines 
the process of simultaneous emission (absorption) of two plasmons by energetic parton 
for the reason discussed in Section 2. The last line of equation ()7.13|1 is a consequence of 
a property of probability symmetry (j3.17j) over permutations of external soft momenta 

^^2(v| p,pi) = W2(v| Pi,p). (7.14) 



Now we turn to analysis of the scattering probability u'2(v| p, — pi). This expression 

can be considerably simplify if it will be used a fact that on plasmon mass-shell the 

(1) 

partial coefficient function (v|p, —pi) satisfies an equality ()6.2|1 . Taking into account 
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this property and a structure of propagator ()7.4|) . a scattering probability can be written 
in the form more convenient for subsequent analysis 




T^(v|p,Pi) + Tll(v|p,pi) + T^(v|p,Pi- ' 



where we extracted all kinematical factors from the scattering amplitudes 

T^(v|p,Pi) = Tll(v|p,pi) = *r°°^(p,-Pi,-g)4f^V''*A'(?) 



5 



(7.15) 



V -pi . - - . - - " qo\ ^ 

T^(v|p,Pi) = 5r°o^(p, -Pi, -q){5^'' - ^y'*A\q). (7.16) 

Here, q = p — pi = {u^ — cjp^, p — pi) is the soft energy and momentum transfers. It is 
convenient to interpret the terms T*^, T^' and T-*- by using a quantum language. The term 
T"-^ is connected with the Compton scattering of the soft modes (plasmons) by energetic 
parton. T" defines the scattering of a quantum oscillation through a longitudinal virtual 
wave with propagator *A'(q'), where a vertex of a three- wave interaction is induced by 
three-gluon HTL-amplitude *r'^^). T-*- defines the scattering of a quantum oscillation off 
energetic parton through a transverse virtual wave with propagator *A*(g). Here, a vertex 
of a three-wave interaction is induced by HTL-correction ST^^^ only, since the contribution 
of the bare three-gluon vertex drops out. 

In the subsequent discussion it can be used a general phylosophy, proposed by Braaten 
and Thoma in ^H] applying to this case. Introduce an arbitrary momentum scale |q*| to 
separate the region of soft^^ momentum transfer |q| ~ gT from ultrasoft one |q| ~ g^T. 
It should be chosen so that g^T ^ |q*| <^ gT, |q*| ~ g^T\n{l/g), that is possible 
in the weak-coupling limit for g —>■ 0. The general analysis of contribution to energy 
loss of three scattering amplitudes each on the right-hand side ()7.15j) shows that a basic 
contribution is determined by the last term with its part, that contains a transverse part 
of gluon propagator *A*(g) in a region of momentum transfer |q| < |q*|. This fact caused 
by two reasons. The first one is connected with existence of the infrared singularity in 
integrand ()7.13j) . that generated by absence of screening in scalar propagator *A*(g) for 
Qo "C |q| < |q*|. Entering the magnetic screening "mass" fi in the transverse scalar 
propagator we have logarithmic enhancement in comparison with other contributions 
both in a region |q| < |q*| and in region |q| > |q*|. The second reason is associated 
with existence of some effective angle singularity in a medium induced vertex function 
6T^'^^{p, —pi, —q) with its convolution with a transverse projector {5^^' — q-'q-''/^^)- These 
two facts will be considered in detail in the next Section. 



^^By the fact that we restrict our consideration to study energy loss defined by scattering on longitudinal 
plasma waves, consideration of hard momentum transfer makes no sense. At the hard momentum scale 



plasmon mode is overdamped j47| . 
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8 Energy loss caused by scattering off colorless plas- 
mons 



In a region of small momentum transfer |q| < |q*| the approximation uj^ = t^p^+q — 
cjp^ + q ■ holds with = duj'p_^/ dpi, and such the resonance condition reads 

- ^pi - V ■ (p - pi) ~ u ■ q = 0, 

where u = v— Vp^ is relative velocity. At this point, it is convenient to change in Eq. ()7.13|) 
the integration variables from p, pi to pi and q. In kinematical factors on the right-hand 
side of fl7.15|) we can set p ~ pi and by this means instead of ()7.13|1 we obtain 

(8.1) 

Here, on the right-hand side, by reason outlined in closing previous Section we leave only 
the scattering amplitude T-'-(v|p, pi) fEq. ()7.16p ). As it will be shown in subsequent 
discussion this amplitude contains no singularities for |u| = and |q| = 0, and therefore 
(5-function in integrand ()8.1jl can be written in the form 

5(u ■ q) = 1 — - — - 5(cosa), 
|u| |q| 

where a is angle between vectors u and q. For integration over momentum transfer q it 
is more naturally to introduce the coordinate system in which axis OZ is aligned with the 
relative velocity u (FigEl) for fixed momentum pi. Then the coordinates of vectors q and 
Pi are equal to q = (|q|, a, ■?/'), pi = (|pi|,ai,0). By <I> we denote the angle between pi 
and q: (pi ■ q) = |pi||q| cos$. The angle $ can be expressed as 

cos <l> = sin a sin ai cos ip + cos a cos ai, (8.2) 

and the integration measure is ciq = c[^d\q\dtlj sin a da. The integral over polar angle a by 
virtue of 5-function determines the value a = 7r/2 in the integrand ()8.H) . Futhermore for 
integration over momentum pi it is convenient going to coordinate system in which axis 
OZ is aligned with velocity v of energetic parton, as it was depicted in FigjTI Then the 
coordinate of vector pi is equal to pi = (|pi|, -i?, 0). The following relations 

|u| cosai = |v| cosi? — |Vp|, 

|u| sinai = |v| siwd 
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are true. The integration measure with respect to "pi"-space is dpi = pld\q\d(j)smi)di). 
Taking into account above-mentioned, the expression ()8.1|) can be presented in the form 



dx J 



2tt /Nra. 



2tt 



27r2 
|q1 



' Z;(Pl 
2uoi 



pi 



.3) 



X 



dcos-d- 



(i|q||q| dcosa5{cosa) T"'"(v|p,pi) 







-1 



Running ahead we note that angle singularity mentioned in closing previous Section is con- 
nected with singularity with respect to angle variable ai. This singularity has dynamical 
origin, since it connected with medium induced three- wave vertex 6T^^\ The consequence 
of this angle singularity will be appearing selected directions in initial isotropic hot gluon 
plasma. It is the directions of primary rescattering of plasmons by hard parton, or di- 
rections over which a main energy loss of parton arise. By virtue of this fact, the basic 
characteristic of interaction process of a high energy incident parton with a hot QCD 
medium, in addition to integral energy loss ()8.3|1 . is the angular dependence of the energy 
loss 

/ dE^^Y 
\ dddx ) 

with respect to velocity direction v. We turn to derivation and analysis of this function. 

Using an explicit expression for HTL-correction 5T^^^ (Eq. (3.30) in work Frenkel and 
Taylor [2]), the initial scattering amplitude T-*- can be presented in an analytic form 



lv|p,Pi, 



(v ■ (n X q)) 



-fq^n"'^(g) + 



2n^ 



IP • qjp 



(p ■ qVi - ^(pi ■ q)] 5r°°°(p, 



ipi • qjPi 



on— shell' 



where n = (p x pi), and vertex correction 5T^^^ is 



5V^'^\p, -pi, -q) = -mj, { cuiM(pi, q) - cuM(p, q) + 



tnqo 



-A) 



with 

M(p,g) 



In 



p-q + J-A{p,q) 



and finally gluon self-energy is 



n""(g) 



m 



D 



|q| 



F{x) 



A{p,q) = p q -{p-q) 



on— shell 



< 0, 



In 



1 + x 



X 



i7i9{l - \x\) 



i.6) 
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where qo = (v • q). In approximation of a small momentum transfer |q| < |q*| and large 
phase velocity of plasmons a;py|pi| ^ 1, the expressions in square braces in Eq. ()8.4|) can 
be set equal 

(p ■ q)^^^ - J^ijpi ■ q) ~ q'^^^ , 2vl^ + (p ■ q)p2 - (pi ■ q)p2 ^ Sq'p?, 
and the medium-induced vertex 5T^^^ reads 



D 



go 

2|q| 



In 



|q| + ?o 



|q| - go 



go 
2|q| 



As we see from Eq. fl8.4p . in this approximation the exact cancellation of the contribution 
with (5r°°'' with gluon self-energy ()8.fij) is arisen, and therefore expression for scattering 
amplitude T-*- has a simple form 



T^(v|pi,q)^3prA*(g) 



V- ((qx pi) X q) 
(q X Pi)^ 



We approximate the transverse scalar propagator *A*(g) by its infrared limit as go ^ 
|q| ^ over surface (v ■ q) = (vj,^ ■ q), namely 



*A*(g) 



nx 



(q2 + /i2)-— m^|v|,J cos$ 



where the angle $ is derived by relation ()8.2|) . Here, we introduce the magnetic "mass" 
/i to screen the infrared singularity. 



After integration over angle a in Eq. ()8.Hj) the expression for amplitude ()8.8|1 reads 

lul cosai 



T^(v|pi,q)^3pf*A*(g) 



1 — sin aicos^il) ' 



.10) 



and in the propagator ()8.9|) it is necessary to set cos$ = sincti cos'?/'. The integrand 
for loss ()8.3|) represents a very complicated function of angle ip. However, here, one can 
use the fact that propagator ()8.9p represents a "smooth" function of angle variables with 
respect to the function going from vertex part in Eq. ()8.1()j) . The last one has a singularity 
with respect to variable cosai. For this reason within accepted computing accuracy we 
can factor out the function |*A*(g)p from the integral with respect to angle ■?/', replacing 
it to average over -0 on interval (0, 27r) 



W(g)p 



|*A*(g)|2 = -/rA*(g)p# 



(q2 + ;u2) J(q2 + /i2) + ^piSin^"! 



f8.111 
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where Xpi = (^/2)^D|VpJ. The remainning integral over dip in a general expression ()8.3|) 
is simple calculated and by virtue of (jS.lOp equals 



2lT 



,2, 



dip 1 + cos^ai 



cos ai / -; r—r = TT 







1 — sin aicos'^ip)'^ \ cosai| 













IpiI 


< 


IpiI 






[ 0, 


IpiI 


> 


IpiI 



We remind that angle ai is connected with angle by means of equality |u| cosai = 
|v| cos'd — |VpJ, and therefore the last relation has singularity in point |VpJ = |v| cos-;?. 
This, in particular, enables us to approximate the relative velocity 

|up = + Vp^^ - 2|v||VpJ cos?? ^ v^sin^-(9. 

We set the plasmon number density A^p^ by step function 

where Nq is a certain constant independing of coupling g, parameter 1 < p < 2. The 
parameter p, as it have been discussed in Section 5, determines a level of soft plasma 
excitations (or value of plasmon occupation number). In deciding on value |pi|max we can 
guide by the next circumstance. The leading order dispersion equation (1.2.1) defines a 
spectrum of longitudinal oscillations for the entire {u, |p|)-plane. However it was shown 
(for QCD case in Refs. |17j and for SQED case in Ref. P3|), that by virtue of specific 
behaviour of longitudinal oscillation spectrum near light-cone u = |p|, consideration 
next-to- leading terms in the dispersion equation leads to important qualitative change. 
There is only a finite range of |p| in which longitudinal plasmons can exist. Therefore 
as |pi|max one can choose a value of plasmon momentum such that modified dispersion 
curve hit the light-cone. For pure gluon plasma we have 



I |2 o 2 

IPllmax = 3W 



max pi 




1-^ +^ 

2 ^ C(2) 



(8.12) 



where 7^; is an Euler's constant and ( is the Riemann zeta function. 



The chosen approximation for plasmon number density is very crude, and we pursue 
only an aim of obtaining explicit analytical expression for energy loss. The dependence 
A^p on p is determined by solution of kinetic equation p.lj) . where the right-hand side 
represents a sum of collision terms ()2.2|) - ()2.4|) . (1.2.3) - (1.2.5) etc. (see, e.g. for ordinary 
plasma Refs. |1S1)- In this case as an alternative value of Ipilmax, Eq- ()8.12|) . one can 
consider 
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i.e. |pi|max will be dependent on temperature by very nontrivial manner. 

With regard to all above-mentioned approximation and Z'(pi) ^- — > 1, we find the next 
expression for the angular dependence of the energy loss of a hard parton from Eq. ()8.Hj) 



' - ■ ' ^ cu^iATn V sin^t? 



d\pi 




1 pi 




V 


cos-i? 



J d\q\ |q| |*A*(g)|2. 



(8.13) 

Using expression for mean value of scalar propagator (jS.llll . the integral over d\q\ is easily 
calculated and in two limiting cases it equals 

|q1 [ ^ lnP^*n 

l/|q.p V^^^ ^Pl ,/i'>Xpi- 

By using this fact, we can approximate this integral by the following expression 

111 , /, *|2\ . 

I rf|q| |q| I *A*(g)|2 ^ 0(xp, - /i^) ^ l^^^f J + ' ^p^) (§.14) 

Here, on the right-hand side we leave the terms which are more singular over /x^. Magnetic 
screening mass in the approximation ()8.14|) carries out separation of plasmons over group 
velocity on "fast" (xpi > /^^) and "slow" ones (xpi < /i^)- 

We going in Eq. ()8.13|1 from integrating over |pi| to integrating over |VpJ, setting for 
small |pi|: d\pi\ ^ (5/3)co'pi(i|VpJ = (5/3)c<jpi(if 1. Then with regard to ()8.14j) we will have 
instead of 



— — ~ — ^ V sm¥ 



dddx 



^.15) 



x<^inp^ / eU--^)— ^ + ^ / d[-^-vA- 

Y yU^ y J ^ 7rm|)^ fi |fi — |v| cos?7| J ^Trmj-, ' |f 1 — |v| cosi/l 

2 Svr^ 1 2 7>T 3|pi|max 

<^ =— -m^No, vim!,^ = — . 

3 gp bujpi 

The integrals on the right-hand side of Eq. ()8.15|1 are well defined in regions |v| cos'i? < 
and Ivjcos-^ > fimax- In the region < Ivjcos-i? < fimax these integrals contain 
non-integrating singularity, and here, the procedure of their extension to finite value is 
required. 
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Now we introduce the notation for integrals entering into the right-hand side of 
Eq. dHISD 



V(z\ a,b)= [ , , , < a < 6 < 1 ; W(z\ a,b) ^ [ , , , < a < 6 < 1, 

{\y-z\ ly\y-'A 

\z\ < 1. 

We extend the functions V{z\ a,b) and W{z\ a,b) based on requirement of continuous of 
total loss {dE^^^^/dxfP as velocity function of energetic parton in the neighborhood of 
value |v| = "Wimax- The result of extension is 

V{z\a,b) = e{a-z) ln(^^ ) + e{z - a)e{b - z) Inf ^ 



a — zl \{z — a)(b — z) 



+ 9{z - b) In 



z 



W{z\a,b) = 9{a 




+ 9{z - 6) - In 



z \aiz — b) J 

These extended expressions are logarithmic diverge on the endpoint of integration interval. 

The existence of singularities over variable |v| cosi) leads to that the angle distribu- 
tion of energy loss {—dE'^^^^/d'ddxY^ has sharp peak in definite directions with respect to 
vector of velocity v, that as was mentioned above, it can be considered as more probably 
directions of plasmons reradiation, after their scattering off energetic parton. If reradiated 
plasmons are sufficiently energetic, i.e. |VpJ ~ fimax then here, one can tell about forma- 
tion of plasmon jets along these directions, which after hadronization can be manifested 
in the typical hadron distribution. 

Let us calculate the contribution to {—dE^^^^ / dddxY^ of the region of fast plasmons 
^^imax — 5 < |VpJ < max, whcre parameter 5 is defined in the interval < 5 < fimax — 
2/i^/7rm|,. From the general expression ()8.15|1 in this case we have 

~ G Ivl sin 

^Imax '^5'^lmax)? 



d-ddx j^^^^ 



V 27r2 



where the function W is defined by Eq. fl8.17|) . In the Fig.|Hl the dependence of this 
expression from cost!} for 6 = 0,2 and 0,05; |v| = 1 and choice Wimax — 3/5, is given. The 
typical peculiarity of this dependence is presence of sufficiently sharp boundary over angle 
i?, for which the loss become essential. This boundary is determined by cos^q = i^imax/lv]. 
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— d'ddx ) from COS 'i? for 6 — 0,2 (solid) and 5 — 0,05 (point) 

under |v| = 1. 

The angle distribution of energy loss from whole region of rescattering plasmons < 
I'^pJ ^ I'^pjmax, is presented in FigEl It is seen that main loss is connected with creation 
of "slow" plasmons |VpJ ~ and lie in region d ~ 7r/2. 

Now we calculate the total energy loss. Integrating the expression ()8.15j) over angle d 
and using ()8.16|) . ()8.17|) we will have 



+ 



/ dzi\ir\'-z')v(z 0,^) (8.1 

-|v| 



27r2 J |v|2 2^^ 



|v| 



2 ' r / ' 



where in the last line we set fimax = 3/5. Using an explicit expression for function V 
(Eq. ()8.16p ). we derive in limit //^ — an asymptotic expression for contribution from 
"slow" plasmons 



dx I ^„ 21 2 \ 27r2 / \ 2u 



c^s \^q2^i^u const + 0(^2) + _ _ _ ^ ^8 



2 
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226.355 




— Mdx ) from cost? for whole region of rescattering plasmons, 
< |VpJ < Ivpjmax under |v| = 1. 

Similar calculations result in asymptotic expression for the second contribution on the 
right-hand side of ()8.18|) 

We see from Eq. ()8.19|) and ()8.20|) that a main contribution to energy loss in the limit 
/x^ — > is connected with a first term in braces on the right-hand side of Eq. ()8.2()|1 . Thus 
leading contribution to parton energy loss, caused by spontaneous scattering off plasmons 
in hot gluon plasma is defined by expression 

; — ^ — ^ /-''a^T^ n n^ ^ . 8.21 

For highly excited plasma, when p — 2, suppression in coupling constant of this process 
of energy loss as compared with polarization loss is not arisen. Moreover here, 
we have two logarithmic enhacements. The first of them is connected with absence of 
sreening in scalar propagator *A*(g), that gives a factor ln(|q*|yyU^). The second one is 
associated with singularity of integrand of Eq. ()8.3|) with respect to angle variable ai in 
point ai = n/2. This singularity sets a singular curve vi — |v| cost!} = in variables space 
(t?, Vi), whose consequence is appearing of logarithmic singularities with respect to angle "i? 
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in points d ~ 7r/2 and d i3q = arccos (|v|/fimax)- They lead to a factor ln^(7rm|,/2yU^) in 
expession for loss ()8.2H) . It can be assumed that this angle singularity is associated with 
the massless basic constituents of a hot gluon plasma, i.e. with massless of hard transverse 
gluons. The inclusion of the asymptotic thermal mass ^ = g'^T'^Nc/Q for hard thermal 
gluons and use of the improved three-gluon HTL amplitude *r^^j^2(P) ~Pij ~P2) |1H| can 
removes this type of singularity. This hypothesis requires separate and more careful 
consideration. 



9 Off-diagonal contribution to energy loss 

In this Section we briefly analyze on a qualitative level a role of off-diagonal term ()7.9|) in 
the process of energy loss of energetic parton. Gailitis and Tsytovich [HE] have pointed to 
the necessity of considering this contribution within usual plasma, and here, we extend 
the conclusions of this work to non-Abelian case. 

Using the explicit expression for effective current JqI"'^{p) ()7.3|) and leaving only longi- 
tudinal parts in the propagator and correlation function we lead to the expression (instead 

of (uni)) 

X Im(*A'(p))Re (v| p, -pi,Pi)p'p\'pj) 6{v-p). 

The explicit form of the coefficient function K'^'J^^^'^^^{ir\p, —pi, —P2) in linear approxi- 
mation in a color charge Q\ of parton is defined by expression ()4.8|1 . With the help of 
()4.8j) we obtain the required expression, where 

(1) (1) 
K"'^^''''^''^{w\p,-pi,Pi) = Tt{T^T^){K'''''^{v\p,-Pi,pi)+ K<'\v\p,pi, -pi)}. (9.2) 

The special feature of the last expression is the condition that in spite of the fact that 
every term in braces on the right-hand side of equation ()9.2j) contains the singularities 
of the l/(t> ■ p), *D^^'{0) etc. types, however in a sum they are reduced. Therefore the 
expression ()9.2j) is well-definite and equals 

^a6feaui.;(^| _p^^p^) ^ NjA[*r'''''HP,Pl, -P, -Pi) *'D"\p)v'' 

- ( + 7^^*r"^^(p, -pi, -p + pi) *V"'ip - p^)v^' (9.3) 
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-*r''^ {p, pi,-p- pi) *v" {p + pi)v^' 



Futhermore we add off-diagonal contribution with expression for polarization loss 
()7.(ij) . After some regrouping of the terms this sum reads 



rfi5;{o)i\P°i 



dx 



+ 



dx 



pol 



Ai + A2. 



Here, on the right-hand side the function Ai is 



Ai = -^ 

V 



1 (Nm., 



where 



27r2 



dp 



Im CA\p))\l + Re (n(^)'(p)*A'(p)) \6{v ■ p), (9.4) 



p^ 



dpi\ 



X 



^iPij ^1 



{p,pi, -p, -pi)p'p\^p'pi' 



pi 



(9.5) 



and the function Ao has the form 



A, 



N^as 
27r2 



X Im(*A'(p))Re 
1 



dpdpi j 
1 



p'- 



-^(E2),^(v-pi) 



(9.6) 



r"^^ (p, -Pi, -P + J9i) *V"\p - pi)t;^' 



(^ -Pi) 



'r''^{p,pi, -p - Pi) *V"' {p + pi)v^' Wp\'6{v ■ p). 



Let us analyze the expression Ai. It can be shown that the entering function Yi^^^\p) 
represents the correction to longitudinal part of gluon self-energy in HTL-approximation 
(i.e. to n'(p)) taking into account change of dielectric properties of hot gluon plasma 
by the action of the processes of nonlinear interaction (longitudinal) plasma excitations 
among themselves. For low excited state of plasma (5 = 1/2 for estimation ()7.7|) ). corre- 
sponding to level of thermal fluctuation, correction 11*^^^' (p) is suppressed by more power 
of g in comparison with I{\p). In the limiting case of a strong field {5 = 0) the func- 
tion n*^^)'(p) (and also corrections of higher power in {Ef)p) is the same order in g as 
gluon self-energy in HTL-approximation, and therefore considering influence of nonlinear 
plasma processes on the energy loss of energetic parton becomes necessary. 

One can define an effective scalar propagator *A\p), which takes into account addi- 
tional contributions, considering nonlinear effects of plasmon self-interactions. Here, we 
have 

00 

*A-^'{p) = *A-^'{p) - ^ H('^)'(p), (9.7) 



n=l 



49 



where 

The algorithm of calculation of the effective subamplitudes *r'^") entering in integrand, 
was defined in [J. For low excited state, when */Sr^\p) ^ 11'^-'^^' (p) ^ Y[^'^^\p) . . . from the 
equation ()9.7j) to the first order in (E^)p. we obtain an approximation 

whose imaginary part is equal to 
Im {*R{p)) - Im (*A^(p))(l + Re {U^'^\pyA\p))] + Re (*A'(p))Im (n(^)'(p)*A'(p)), (9.8) 



correspondingly. Here, the first term on the right-hand side is precisely coincident with 
appropriate function in integrand of ()9.4|) . If the second term is absent on the right-hand 
side of Eq. fl9.8p . than the contribution Ai to the energy loss can be interpretated as 
polarization loss with regard to change of dielectric property of medium induced by plas- 
mon self-nteraction. This can be effective present as a replacement of the HTL resummed 
longitudinal propagator, entering into expression for polarization loss ()7.(ij) by an effective 
one 

*A'(p) ^ *a'(p). 

However the existence of the second term on the right-hand side of Eq. ()9.8p gives no 
way completely restricted ourselves by such simple transition. The physical meaning of 
this fact is not clear. 

Let us analyzed the contribution A2 (Eq. ()9.fi|l ). This contribution (as previous one) 
is not equil to zero only for plasma excitations lying off mass-shell. From the expression 
for basic diagonal contribution to energy loss ()7.1H) . ()7.12p we see that it contains the 
singularities off-plasmon shell in the forms l/{v ■ piY and l/{v ■ pi), when frequency and 
momentum of plasma excitations approach to " Cherenkov cone" 

{v -pi) 0. 

The existence of these singularities causes the large rotation velocity of color vector of 
energetic parton in random plasma field. However these singularities are exactly compen- 
sated by analogous ones in expressions for energy loss A2, Eq. ()9.(i|l . It is not difficult to 
check this by correlating the expressions ()7.11|) and ()9.6|) . Thus the complete sum of all 
contributions to energy loss of energetic parton for scattering off plasma excitations lying 
off-plasmon shell is a finite value, but as supposed it is sufficiently small in comparison 
with energy loss for scattering off plasma excitations lying on-plasmon shell to neglect it. 
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10 The Fokker-Planck equation for beam of energetic 
partons 

The goal of this Section is a derivation of the Fokker-Planck equation, describing an 
evolution of a distribution function f(k,x,Q) for a beam of energetic color partons^^ 
due to the processes of energy loss and diffusion in momentum space for scattering off 
soft longitudinal plasma excitations. Let us supposed that particles number density is so 
small that one can neglect by partons interaction in beam among themselves. One of the 
consequences of derivation of this equation having relevance to the subject of our study 
will be obtaining more complete expression for energy loss and also a possible gain of 
energy of isolated energetic parton of beam which follows from well-known connection of 
the drag coefficient A*(k) with expression for energy loss dE/dx [HZj, |HHj . 

First of all we consider a more general problem for initial formulation, which in prin- 
ciple gives a possibility to consider a diffusion of beam of energetic color particles in an 
effective color space on an equal footing with diffusion in momentum space. We assume 
that a time-space dependent external color current j^^*"(a;) starts acting on the system. 
In this case the soft-gauge field develops an expectation value = A^^^i^x) ^ 0, 

and the number density of the plasmons N^^^ acquires a nondiagonal color structure. The 
equation describing the time-space evolution of the distribution function /(k, x, Q) in first 
approximation in a number of emitted and absorbed plasmons reads 

I + V ■ 1^ + ^^Q-^ {£\x) + (v X B\x))) ■ /(k, X, Q) 

+ <7r'"^^M«(x)^(Q7(k,x,Q)) (10.1) 

1 r dp dpi {a'},{a}/. /I N| X 

x{Nl^'^iNl^ + ir<fik,x,Q) - (iV;, + l)-'iVK'^V(k + q,x,Q)} 

1 f dp dpi {a},{a'}/, /, x, X 

-27(2^(2^"^ ^(v(k)|pi,p) 

y<{NlfiNl + i)<«^/(k,x,g) - (Nl^ + ir'«iv^^i</(k - q,x,g)}. 

Here, we use a multiindex notation {a} = (a,ai). In deriving collision integrals on the 
right-hand side of Eq. ()1().1|1 we take into account a property of scattering probability 
relative to permutation of the external soft momenta p and pi: 

^j"'^'^"^(v(k)|p , pi) = wt^'^'''\ir{k + q)| pi, p). (10.2) 

^^We will denote the momentum of energetic partons of beam by the same symbol k as momentum of 
hard thermal particles of hot medium (unlike velocity) , lest not introduce a new notation. 
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This property of the scattering probabihty is a generahzation of the property ()7.14|) . 
since the last one is vahd only in the framework of HTL-approximation and for condition 
|q| ^ |k|. The equality ()10.2|1 is a consequence of detailed balancing principle and 
therefore is exact. The existence of two collision terms on the right-hand side of Eq. (jlO.ll) 
takes into account two possible processes of exchange of the beam state of energetic color 
particles. 

In the general case the scattering probability in integrand on the right-hand side of 
Eq. also depends on a color vector Q = (Q^). Suppose that this dependence is 

analytic, i.e. the function wj^" ^'^"^ can be represented as an expansion in powers of Q and 
a first term in this expansion can be written in the form 

'^■^'^^(v(k)l p, pi) = (T-^^(v(k)l p, -pi)Q^) (T*»'<^'(v(k)| p, -pOQ"' 



x27i6{E^-E^+^ + ujI-ujIJ. (10.3) 

Here, E]^ is energy of hard parton and T"'^^''(v(k)| p, — pi) = — z/'*"^^T(v(k)| p, — pi) is 
interaction matrix element. The scalar function T(v(k)| p, — pi) is reduced to the function 
()3.18p in HTL-approximation and small momentum transfers. 

Furthermore we restricted our consideration to study of the scattering process of a 
bundle hard partons off colorless part of soft-gluon excitations, i.e. set on the right-hand 
side of Eq. fITTnil N^^^' = 6"''' N^. If one neglect by influence of a mean field (setting it a 
small), then kinetic equation (10.1) can be written as equation for momentum of zeroth 
order in color charge of distribution function 



/(k,x)= JdQfik,x,Q). 



With regard to a structure of scattering probability ()10.3p . after some algebraic transfor- 
mation and regrouping the terms, initial equation po.lj) can be introduced in the form 
more convenient for subsequent study 




where 

(dfY_ fN^aA'fC 



-vr I 11^1 I dpdpi{w2iir(k)\p,pi)[Nlf{k,x) - NlJ(k + q,x, 

(10.5) 

+ W2mk) I pi, p) {NlJ{k, x) - Nlf{k - q, x)) }5iE^ - i^k+q + uj^ 
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is collision term linear in plasmon number density A^p corresponding to spontaneous scat- 
tering processes^^, and 



TT 



dt V 277 



^) j dpdpi{t/;2(v(k) I p, pi) (/(k, x) - /(k + q, x)) (10.6) 



+ w;2(v(k)| pi, p) (/(k, x) - /(k - q, x))]n^N^,5{E^ - EkH 



is collision term quadratic in N^^ corresponding to stimulated scattering process. Here, a 
constant C2 = Q'^Q"' is the second order Casimir. We have C2 = Ca for beam of high 
energy gluons and C2 = Cp for quarks. In the scattering probability W2(v(k)| p, pi) = 
|T(v(k)| p, — pi)p we take out a coupling constant and isotopic factors in a coefficient 
before integral, and the 5-function expressing the energy conservation - in integrand. 

Let us consider approximation, when exchange of plasmon momentum in scattering 
process is much smaller than momentum of energetic parton from beam, that is |q| = 
|p — Pi| ^ |k|. It enables us to expand the distribution function in powers of |q|: 

/(k ± q, x) ^ /(k, x) ± (^q ■ /(k, ^) + ^ (q ■ 1^) /(k, x)- (10.7) 

The expansion of scattering probability 11)2 requires a larger accuracy. Let us consider 
equality W2(v(k)| pi, p) = W2(v(k — q)|p,Pi), which is a consequence of p0.2|) . We 
expand the right-hand side in a small parameter |q|/|k|: 

w'2(v(k)| pi,p) = w;2(v(k)|p,pi) - ^q- J-^^i;2(v(k)|p,pi) + ... . (10.8) 

The expression t(;2(v(k)| p, pi) is also required to expand in a small parameter |q|/|k|: 

t/;2(v(k) I p, pi) = w!^^ (v(k) I p, pi) + wi''^ (v(k) I p, pi) + . . . . (10.9) 

Substituting the last expansion in the either side of Eq. (jlO.Hj) . we obtain 

t.i°)(v(k)|p,pi) = ..i°)(v(k)|pi,p), (10.10) 

(i)r^nM ^ (i)^-nM A, ^^2°^(v(k)| p,pi) /mm 
W2 (v(k)|p,pi) = ^2' ^(v(k)|pi,p) + q , (10.11) 

etc. The principle term of expansion t(;2°'*(v(k)|p, pi) represents the scattering probability 
in quasiclassical approximation with neglecting of recoil of hard parton. Therefore this 
function is identified with a scattering probability in HTL-approximation (Eqs. ()7.15p . 



16 



Here, we follows by terminology accepted in |50|. 
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()7.16|) ). having the same property of symmetry relative to permutation of soft external 
momenta (Eq. (17^1 ) as 

At first we consider expansion of collision term defining a process of stimulated scatter- 
ing. Substituting an expansion ()l().7j) . into and using the properties 
and (jlO.lH) . we obtain that {df/dty^ can be introduced in the form of a Fokker-Planck 
containing the scattering probability in HTL-approximation only 

where the drag and diffusion coefficients respectively are 

A^*'^(k) = -TT (^)'(^)/g^(q- l^^f (v(k)|p,pO)iVpiVp,5(t;- (p-pi))rfpd^^ 

(10.12) 

Bst,^i(k) = n I (v(k)l p, p,)N^N^, 5{v ■ (p - p^))dpdp,. 

In the case of collision term, associated with spontaneous scattering process (Eq. ()10.5|l ) 
a situation is more complicated. If we perform calculations similar to previous case then 
we will have 

where the drag and diffusion coefficient are 

Bsp,.,^j^) = I (^)'(^)/?V^f (v(k)|p,pi)(iVi + iViJ5(^ . (p-pi))rfprfpi, 
and the function $*-^^[/] is defined by expression 

^'''[f] = I (^)'(f )/{-i^Hv(k)|p,pO + ^i^^(v(k)|p,,p)}(Ari- AT^^ 

X • 5{y-{p- pi))dpdpi. 

The existence of the last term $''^''[/] containing a correction to quasiclassical scattering 
probability shows that for collision term connected with spontaneous scattering process 
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in limit of a small momentum transfer and accepted calculation accuracy, it is impossible 
to close on quasiclassical description level. In other words, it is not sufficient to use the 
properties (jlO.f 0|) and ()10.11|1 to remove by the next term in expansion of the scattering 
probability ()1().9|) considering in fact an effect of quantum response of hard parton. 

Nevetheless, having an explicit expression for coefficient functions ()10.12|) . ()10.13|) we 
can write the expression for full drag coefficient 

A'(k) = A^P'^(k) + A'^''^(k) 

and in so doing to define the energy loss by using relation jHZj, jHH] 

The first term on the right-hand side exactly reproduces the expression for energy loss 
()7.13|) used in previous Sections for C2 = Ca- The second term defines acceleration 
obtained by hard color particle for Compton plasmon scattering. As it is seen from this 
expression in fact it equals the product of total section of this scattering process in HTL 
approximation by mean energy of longitudinal QGP excitations. It is necessary to take 
into account this contribution in complete balance of energy loss. The third term is 
formally suppressed in comparison with the first factor |q|/|k| -C 1. However it can be 
expected, that for highly excited state of system, when the plasmon number density is very 
large, the term in integrand proportional to iVpiVp^ (stimulated scattering process), can 
give appreciable contribution to general expression for energy loss. By virtue of the fact 
that this contribution is indefinit in a general case then here, both drag and acceleration of 
partons from beam are possible. In particular it depends on concrete form of the plasmon 
number density A^p as a function of soft momentum p. 

11 Conclusion 

In conclusion we would like to dwell on some important moments which remain to be 
touched upon in this work and on the whole concerned with a mechanism of energy loss 



55 



studied in previous sections. 

In studies of plasmon-hard particle scattering processes and energy losses their related, 
we consider that a mean-gluon field .4^ (a:) in a system equils zero or is so weak that it can 
be neglected by its infiuence on dynamics of the processes in QGP^^. This in particular 
leads to the fact that we have considered interaction of energetic color particle (or thermal 
test particle) with colorless part of soft-gluon excitations only. However, opposite case 
- very strong mean gluon field that is to be necessary account exactly in all orders of 
perturbation theory, is of great interest. In this case, for example, coefficient functions 
p, — pi, . . . , — Ps) under the integral sign of effective currents ()7.Hj) depend by 
themselves on ^^(x), i.e. they can contain insertions of an external or self-consisting 
mean field, in general, of arbitrary order. The interest to such setting up a problem is 
specified by that it is simpler model for research of jet quenching phenomenon in heavy ion 
collisions, whereby high-pt partons get depleted through strong (classical) color field. The 
latter is encountered in the Color Glass Condensate (CGC) jnij. About ten years ago 
Mc.Lerran and Venugopalan [22] suggest that the prompt phase of heavy ion collisions 
is CGC, and therefore infiuence of CGC on energy losses of color partons can be very 
essential. The example of considering CGC in calculation of energy loss has been given 
in Introduction. This is a mechanism of synchrotron-like radiation of an ultrarelativistic 
charge in QCD going through a constant chromomagnetic field, proposed by Shuryak and 
Zahed j2Dl- 

Another important moment that was not touched in this work is concerned with 
consideration of contributions to energy losses of terms higher in powers of the field in 
expansion of the effective current ()7.2|1 . s > 3. In principal, in limiting case of a strong 
random field |A^°^°(p)| ~ l/g{gTY fEq. ()7.7|) ) account must be taken of the whole of 
series in expansion of the effective current, since formally each term of the series becomes 
of the same order as a first term JqJ^. A preliminary analysis of contribution to energy 
loss of higher scattering processes shows that in each subsequent 'diagonal' contribution 
(—dE^^'ydxY'^, n > 2, a term having practically the same peculiarity as {—dE^^^/dxY^ 
can be separated. However, here, we come up against the further problem not allowing 
directly summarized all relevant terms. These contributions contain the pathological 
terms generated by product of two 5-functions of [5(Re *A~"(p + pi))]^ type etc. This 
strongly remind a situation with that one runs in out of equilibrium thermal field theory 
connected with so-called pinch singularity [SSj, with the distinction that argument of 5- 
function contains inverse resumming propagator instead of bare one. For this reason a 
problem of contribution to energy loss of higher spontaneous scattering processes requires 

^^The exception here, is only the beginning of Section 10. The kinetic equation 110. 1|1 written there 
represents a simplest example of that as far as a dynamics of a system for A'^{x) ^ becomes more 
complicated and intricated. 
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separate and careful study. 

A futher remark is concerned with existence of soft-gluon transverse excitations in the 
system. Entering the transverse oscillations in consideration result in appearing additional 
channels of energy losses of energetic particle similar to ones that we study for longitudinal 
oscillations. Besides the scattering processes off pure transverse plasma waves here, it 
is also possible the scattering processes resulting in change of polarization type of soft 
excitations (so-called a conversion processes of plasmons to transverse waves^^) and more 
complicated processes being a combination of usual scattering process and the fusion 
process of two plasmons into one transverse eigenwave, as depicted in Fig.El These 




Figure 10: A scattering-fusion process of two soft-gluon longitudinal plasma excitations off 
energetic color particle with transition to one soft-gluon transverse plasma excitation. 

processes also need to be taken into account for total balance of energy losses. However, 
here, a new specific effect, namely for transverse mode is appeared. This is a possibility of 
large energy-momentum transfer from ultrarelativistic particle for scattering off plasmons 
with its conversion to transverse quantum of oscillations (see the last footnote). The 
energy of secondary quantum may be close to energy of initial hard particle. This is 
connected with that unlike longitudinal excitations, transverse ones have no restriction to 
maximum possible value of momentum of ()8.12j) type. However when energy of secondary 
quantum is close to energy of fast particle, then quantum effects become essential. Here, 
it is necessary to use quantum expressions for scattering processes on gluons of Compton's 
type with change of incoming line of vacuum gluon by plasma gluon with longitudinal 
polarization, and a factor provides renormalization of a gluon wave function by thermal 
effects. The scattering on longitudinal gluons with conversion to hard transverse radiation 
produces additional mechanism imposed upon gluon bremsstrahlung. The processes of 
such a kind for usual (Abelian) plasma were first considered by Gailitis and Tsytovich in 
Ref. . 

^*For example, on Fig.^and Fig. |21 incoming soft-gluon line have a longitudinal polarization, whereas 
outgoing soft gluon line - transverse ones. 
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Finally, a last remark that should be mentioned, is concerned with energy loss of 
color dipole, next in order of complexity to color object after single color-charge particle. 
Classical color dipole can be considered as 'toy' model for heavy quarkonia (cc or bb) 
propagating through QGP, and its energy loss may be relevant for the distortion of heavy- 
quarkonium production by plasma formation. Within linear response theory this problem 
was first considered by Koike and Matsui in Ref. ^Tj. However an attempt of direct 
extension of the result obtained by Koike and Matsui to the case of the scattering of color 
dipole off soft-gluon excitations of QGP encounter specific difficulties. In fact we can 
easily define expression for "starting" current generated by the color pairs. For a neutral 
pair, in the approximation 'point' dipole we have 

J^tip) = (^^"^(^ ■ P) ^{^^^^ (P ■ ■ v) + ^ (P X V)(d" X v)}, (11.1) 

where d** = (Q"q — Q2o)^o is the color dipole moment with fixed separation ro of color 
charges in a intrinsic frame of reference, and v is a velocity of dipole. Substituting pi.l|) 
into expression for energy loss ()7.1|1 . previously replacing 

JdQ^J dQ.dQ,, E^\p) ^ E-q(p) = -too *V^^{p)J^qM^'^]{p), 

we derive polarization energy loss for color dipole. To take into account scattering process 
for color dipole off soft plasma waves and energy loss associated with this process, we must 
calculate the next term in the expansion of effective color-dipole current JgQ[A^^^{p), 
linear over soft fluctuation field of system A^^^"". The problem lies in the fact that not only 
a fleld of incident wave, but also field of second charge of pair are acted on each of color 
charge vectors of dipole. This circumstance very complicates and tangles the problem, 
compared to single color charge (see, e.g. work of Khriplovich |55| on construction of a field 
produced by two fixed Yang-Mills point charges), and therefore demands development of 
more refined method for calculating effective color-dipole current, that is out of the scope 
of the present study. 
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